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ABSTRACT 

Aims. To measure the supernova (SN) rates at intermediate redshift we performed the Southern inTermediate Redshift ESO Supernova 
Search (STRESS). Unlike most of the current high redshift SN searches, this survey was specifically designed to estimate the rate for 
both type la and core collapse (CC) SNe. 

Methods. We counted the SNe discovered in a selected galaxy sample measuring SN rate per unit blue band luminosity. Our analysis 
is based on a sample of ~ 43000 galaxies and on 25 spectroscopically confirmed SNe plus 64 selected SN candidates. Our approach is 
aimed at obtaining a direct comparison of the high redshift and local rates and at investigating the dependence of the rates on specific 
galaxy properties, most notably their colour. 

Results. The type la SN rate, at mean redshift z = 0.3, amounts to 0.22^°+° •}* h 2 10 SNu, while the CC SN rate, at z = 0.21, is 
0.82^24^ 26 ^70 SNu. The quoted errors are the statistical and systematic uncertainties. 

Conclusions. With respect to local value, the CC SN rate at z = 0.2 is higher by a factor of ~ 2 already at redshift , whereas the type 
la SN rate remains almost constant. This implies that a significant fraction of SN la progenitors has a lifetime longer than 2-3 Gyr. 
We also measured the SN rates in the red and blue galaxies and found that the SN la rate seems to be constant in galaxies of different 
colour, whereas the CC SN rate seems to peak in blue galaxies, as in the local Universe. 

SN rates per unit volume were found to be consistent with other measurements showing a steeper evolution with redshift for CC SNe 
with respect to SNe la. 

Finally we have exploited the link between SFH and SN rates to predict the evolutionary behaviour of the SN rates and compare 
it with the path indicated by observations. We conclude that in order to constrain the mass range of CC SN progenitors and SN la 
progenitor models it is necessary to reduce the uncertainties in the cosmic SFH. In addition it is important to apply a consistent dust 
extinction correction both to SF and to CC SN rate and to measure SN la rate in star forming and in passive evolving galaxies in a 
wide redshift range. 

Key words. supernovae:general - star: formation - galaxy:evolution - galaxy:stellar content 



1. Introduction 

SNe, the catastrophic explosions that terminate the life of some 
stars, play a pivotal role in several astrophysical topics. They 
provide a crucial test for stellar evolution theory and, as the main 
contributors of heavy elements in the Universe, they are key ac- 
tors of the chemical enrichment of galaxies. Their rapidly ex- 

Send offprint requests to: M.T. Botticella, e-mail: 
botticella@oa-teramo . inaf . it 

* Based on observations collected at the European Southern 
Observatory, using the 2.2m MPG/ES O telescope on the La Silla 
(ESO Programmes 62.H-0833, 63.H-0322, 64.H-0390, 67.D-0422, 
68.D-0273, 69.D-0453, 72.D-0670, 72.D-0745, 73.D-0670, 74.A-9008, 
75.D-0662) and using Very Large Telescope on the Cerro Paranal (ESO 
Programme 74.D-0714). 



panding ejecta sweep, compress and heat the interstellar medium 
causing gas outflows from galaxies and triggering the star forma- 
tion process. Furthermore, in the last few years type la SNe have 
consolidated their prominence as cosmological probes provid- 
ing the first direct evidence for the acce leration of the Universe 
dPerlmutter et al.lll999HRiess et alj|1998l) . 

SN statistics is another important cosmological probe even 
though less exploited. In particular the SN rate as a function of 
the cosmic time is linked to some of the basic ingredients of the 
galaxy evolution such as mass, SFH, metallicity and environ- 
ment. 

The rate of CC SNe, including type II and type Ib/c, gives, 
for an assumed initial mass function (IMF), a direct measure- 
ment of the on-going SFR because of the short evolutionary life- 
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times of their progenitors. On the other hand the rate of SNe la 
echoes the long term SFH because these SNe, originating from 
low mass stars in binary systems, are characterized by a wide 
rang e of delay times between progenitor form ation and explo- 
sion ( Yungelso n& Livio Il2000t iGreggio 112005b . 

Conversely, if a SFH is assumed, the progenitor sce- 
narios for SN la can be constrained measuri ng their rate 
as a fu nction of galaxy ty pe and re dshift dMadau et al. I 
1998t iDahlen & F ransson Il999t iPainetall 120021: 

iTonrv et al. I 
2005t 



Strolger et alTT2004t 
2003t IScannaoieco & Bildsten 



Maoz & Gal-Yam 

"""20051 iMannucci etal 



2004; 



Man nucci et al. I I2006L iNeill et al. 1 120061: ISullivan et al71 
Forster et al] (2006J). Despite promising prospects, measure- 
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ments of SN rates were very scanty so far. 

A new opportunity for SN searches emerged in recent years 
thanks to the availability of panoramic CCD mosaic cameras 
mounted at medium/large size telescopes. These instruments al- 
low the deep monitoring of wide areas in the sky and thus the 
collection of statistically significant samples of SNe at interme- 
diate and high redshifts. 

In this paper we report the results of the "Southern inTerme- 
diate Redshift ESO Supernova Search" (STRESS) which was 
carried out with ESO telescopes and designed to search for SNe 
in the range 0.05 < z < 0.6. 

The paper is organized as follows: in Sect. 2 we illustrate the 
aims and the strategy of STRESS, in Sect. 3 we present the data 
set collected during our observing programme. The following 
three sections describe the data analysis, namely the selection of 
the galaxy sample (Sect. 4), the detection and classification of 
the SN candidates (Sect. 5) and the estimate of the control time 
of the search (Sect. 6). The measurement of the SN rates and 
the estimate of the statistical and systematic uncertainties are 
presented in Sect. 7. Our results are compared with published 
measurements and models predictions in Sect. 8. Finally, a brief 
summary is given in Sect. 9. 

We remark that the measurements reporte d in this paper su- 
persed e the preliminary results published in Capp ellaro et al. I 
(2005) because of i) better SN statistics , if) better filter cov- 
erage for characterization of the galaxy sample, Hi) improved 
correction for the host galaxy extinction. 

Thorough the paper we adopt the cosmological parameters 
H<£\ = 70kms~ 1 Mpc~ 1 , Q M = 0.3, Q A = 0.7. Magnitudes are in 
the Vega System. 



2. Search goals and strategy 

Most past and current high-redshift SN searches were/are aimed 
at detecting SNe la for measuring cosmological parameters, 
while the measurement of SN rates is a secondary goal. As such, 
the strateg y of these searches are affected by seve r al shortcom- 
ings (eg. ISchm idt et al.~l ll998t iPain et al. 1 120021: ITonrv et al. I 
2003: lDahlen et al. Il2004t INeill et al. 112006b . In particular, high- 
redshift SN searches are optimized for the detection of un- 
reddened SNe la near maximum light, so that all variable sources 
unlikely to fall into this category are assigned a lower priority 
for follow-up, when not completely neglected. This strategy in- 
troduces a severe bias in the SN sample; besides, in general, it 
prevents a reliable estimate of the rate of CC SNe. 

STRESS was specifically designed to measure the rate of 
both SNe la and CC SNe at intermediate redshift, to be compared 
with the local values and estimate their evolution. To this aim we 



tried to reduce as much as possible the biases which affect the 
different SN types. 

To measure the SN rate we counted the events observed in a 
selected galaxy sample. This is the same approach used to derive 
the local rates, and offers two advantages. By preserving the link 
between the observed SNe and the monitored galaxies we can 
investigate the dependence of the rates on specific galaxy prop- 
erties, most notably their colours, which depend primarily on the 
galaxy SFH, and on metallicity and dust extinction. In addition, 
the SN rates measured per unit of galaxy B luminosity, that is in 
SN unit (SNvfl), can be directly compared with the local ones. 

For the selection and the characterization of the galaxy sam- 
ple, we used multi-band observations and the photometric red- 
shift technique (photo-z), which allow us to derive the redshift, 
the absolute luminosity, and the rest frame colours for all the 
galaxies in our sample. 

In summary our approach involves the following steps: 

1 . the selection of the galaxy sample and its characterization, 

2. the detection and classification of SN candidates, 

3. the estimate of the control time of the galaxy sample, 

4. the measurement of SN rates, the analysis of their depen- 
dence on the galaxy colours, and their evolution with red- 
shift. 



3. The observing programme 

Observations were carried out using the Wide Field Imager 
(WFI) at the 2.2m MPG/ESO telescope at La Silla, Chile. WFI is 
a mosaic camera consisting of 2 x 4 CCDs, each of 2048 x 4096 
pixels, with a pixel scale of 0.238 arcsec and a field of view of 
34 X 33 arcmin 2 . The individual chips are separated by gaps of 
23.8 arcsec and 14.3 arcsec along right ascension and declination 
respectively, for a resulting filling factor of 95.9%. 

We performed observations in the B,V,R,I bands using 
the following ESO/WFI broad-band filters: B/99, fi/123, 7/89, 
tfc/162, Ic/lwp. 

Our initial list included 21 sky-fields evenly distributed in 
right ascension to reduce the observing scheduling requirements. 
The fields were chosen at high galactic latitude to reduce stellar 
crowding and galactic extinction and with few bright sources to 
minimize CCD saturation and ghost effects. 

The observing programme, distributed over a period of 6 
years, from 1999 to 2005, can be divided into three phases. 
During the first year we carried out a pilot programme aimed 
at tuning the observing strategy and testing our software (at the 
beginning of our programme not all filters were available). In a 
second phase (2001-2003, 12 observing runs) we performed the 
SN search in the V band targeting events at redshift z ~ 0.2-0.3. 
In the last phase (2004-2005, 4 observing runs), we performed 
the search in R band aiming to detect SNe at higher redshifts. A 
total of 34 nights were allocated to our programme out of which 
9 were lost due to bad weather. Seeing was in the range 0.7 -1.5 
arcsec, with an average around 1 arcsec. 

The temporal sampling, on average one observation every 3- 
4 months, was tuned to maximize event statistics. The observing 
log is reported in Tab [2] 

Typically, in each run we observed the same fields for two 
consecutive nights, in the first one with the search filter, V or R, 
and in the second one with a different filter to obtain colour infor- 
mation both for SN candidates and galaxies. Each observation, 
typically of 45 min integration time, was split in a sequence of 



#o/70kms-'Mpcr 



SNu = 5A'(100yr)- 1 (lO 10 ^)- 
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Table 1. List of the fields with multi-band coverage. 



4. The galaxy sample 



Name 


R.A. (J2000.0) 


Dec. (J2000.0) 


E(B-V) # 


Band 


03Z3 


03:39:31.1 


-00:07:13 


0.081 


B,V,R,l 


10Z2 


10:46:45.8 


-00:10:03 


0.040 


B,V,R 


13Z3 


13:44:28.3 


-00:07:47 


0.026 


B,V,R 


22Z3 


22:05:05.8 


-18:34:33 


0.026 


B,V,R 


AXAF 


03:32:23.7 


-27:55:52 


0.008 


B,V,R,I 


EF0427 


04:29:10.7 


-36:18:11 


0.022 


B,V,R,I 


EisAl 


22:43:25.6 


-40:09:11 


0.014 


B,V,R 


EisB 


00:45:25.5 


-29:36:47 


0.018 


B,V,R 


EisC 


05:38:19.8 


-23:50:11 


0.030 


B,V,R,I 


EisD 


09:51:31.4 


-20:59:25 


0.040 


B,V,R,I 


Field2 


19:12:51.9 


-64:16:31 


0.037 


B,V,R 


J1888 


00:57:35.4 


-27:39:16 


0.021 


B,V,R 


white23 


13:52:56.4 


-11:37:03 


0.067 


B,V,R,I 


white27 


14:10:59.7 


-11:47:35 


0.057 


B,V,R,I 


white3 1 


14:20:15.3 


-12:35:43 


0.096 


B,V,R,I 


whitehz2 


13:54:06.1 


-12:30:03 


0.073 


B,V,RJ 



# Galactic extinction from Schlegel et al. ( 1998) 



three exposures jittered by 5-10 arcsec to allow a better removal 
of cosmetic defects, cosmic rays, satellite tracks and fast mov- 
ing objects. Due to technical and weather limitations and, occa- 
sionally, to scheduling constraints, in many cases we could not 
maintain our observational strategy. As a consequence, for many 
SN candidates we are missing colour information and for a few 
fields we have insufficient filter coverage to apply the photo-z 
technique to galaxy sample. 

In order to secure the classification of the SN candidates 
spectroscopic observations were scheduled about one week 
after the search run. In the first year, spectra were taken 
at the ES03.6m telescope equipped with ESO Faint Object 
Spectrograph and Camera (EFOSC); in the following years, 
we used the FOcal Reducer and low dispersion Spectrograph 
(FORS1/2) at the VLT. In total 2.5 nights were allocated at the 
ES03.6m telescope and 12 nights at the VLT, 4 of which were 
lost due to bad weather conditions. 

For a better subtraction of the night sky emission lines we 
selected grisms of moderate resolution (FORS1/2 grism 300V 
and/or 300/) which allow the sampling of a wide wavelength 
range (445-865 nm and 600-1 100 nm respectively) with a reso- 
lution of ~ 10A. The width of the slit was chosen to match the 
seeing but, to reduce the contamination of the host galaxy back- 
ground, never exceeded 1 arcsec. Depending on the magnitude 
of the SN candidate, the exposure time of spectroscopic obser- 
vations ranges from 900 seconds to 3 hours. 

As we will detail in Sect. 5 we could obtain direct spectro- 
scopic observation only for a fraction (~ 40%) of the SN candi- 
dates. In order to obtain the redshifts of the host galaxies of the 
remaining SN candidates, and to check for signs of the presence 
of Active Galactic Nuclei (AGN), we carried out a follow-up 
program using FORS2. The observing strategy, tuned to make 
the best use of non-optimal seeing and sky transparency condi- 
tions, proved to be very successful: we obtained 44 spectra of 
SN candidate host galaxies in a single ESO period. 

In addition, a dedicated service program using WFI (12 
hours) was designed to complete the photometric band coverage 
of the fields monitored during SN search. 

Eventually we were able to secure B, V,R,I imaging for 11 
fields and B,V,R for 5 fields, whose coordinates, galactic ex- 
tinction and band coverage are shown in Table Q] Further data 
analysis has been restricted to these 16 fields. 



In this section we describe the steps required to select our galaxy 
sample, namely the detection of galaxies in the monitored fields 
(Sect. l4~TT i. the selection criteria and the application of the photo- 
z technique (Sect. l4~3T >. We also present a number of crosschecks 
with other observations and galaxy catalogues that were per- 
formed to validate our catalogues (Sect. l4~2l and Sect. 14.3b . 

4.1. Source catalogues 

For optimal source detection in each field we obtained deep im- 
ages in each band. First we selected all the images with the best 
seeing and sky transparency for each field and band. Then we 
selected, for each band, a seeing range of < 0.15 arcsec with 
the requirement that it contains the maximum number of images 
which hereafter are co-added. 

All images were preliminarily processed using the /RA/f3 
package MSCRED, spe cifically designed to handle mosaic CCD 
images dValdeslfl998h . We followed the standard recipes for 
CCD data reduction applying bias subtraction, flat field correc- 
tion, trimming and, for / band images, de-fringing. Astrometric 
calibration for each image was perform ed using as reference the 
USNO-A2 catalog dMonet et all 19971) with the task msccmatch. 
The r.m.s. dispersion of the absolute astrometry is ~ 0.3 arcsec. 

The photometric cal ibration was o btained with reference to 
Landolt standard fields (lLandoltll992l) observed on photometric 
nights. Images obtained under non-photometric conditions were 
scaled to the flux of calibrated images by matching the photom- 
etry of the common stars. 

After a proper match of the flux scale, the selected images of 
a given fi eld and band were stacked together using the SWARR^ 
software (Berti n et al. 1 12002) with a third-order Lanczos kernel 
for the resampling, and a median algorithm for the co-adding 
procedure. 

The limiting magnitudes of the co-added images reflect the 
very different total exposure times for the search and comple- 
mentary bands: the 3-cr limit for point sources was reached at 

B Um ~ 23 mag, Vu m ~ 24 mag, R tim ~ 23 mag and l Um ~ 22 
mag. 

In order to consistently measure source colours, we have to 
take into account that the co-added images of a given field show 
different point spread functions (PSF) in each band. To avoid 
colour offsets, the co-added images in each band were convolved 
with an appropriate Gaussian kernel, using the IRAF task gauss, 
to match the PSF of the co-added image with the worst seeing. 

For source detection and photometry we used SExtractor 
dBertin & Arnoutsll 19961) in dual image mode: source detection 
and classification was obtained from the original V band im- 
ages, since this is the band with the best S /N ratio, while the 
magnitudes of the sources were measured on the convolved im- 
ages using the same adaptive aperture as in the original V im- 
ages. This procedure assures that the fluxes are measured in the 
same physical region of each source in all bands. We adopted a 
source detection threshold of 3<x above the background noise and 
a minimum of 12 connected pixels above the detection threshold. 
To measure magnitude we chose a Rron-like elliptical aperture 
(with 2.5 pixel Kron radius) adaptively scaled to the object di- 
mensions (magjauto). 



IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
4 www.terapix.iap.fr 
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Fig. 1. The stellarity index of the sources detected in one of our 
co-added image as function of magnitude. The vertical line indi- 
cates the cut-off magnitude for our galaxy sample. 



Fig. 2. Differences in the equatorial coordinates between our 
source list for the AXAF field and the COMBO-17 catalog for 
the objects in common. 



The instrumental magnitudes were normalized to 1 sec ex- 
posure time and corrected for the instrument zero point, atmo- 
spheric and Galactic extinction. For our purposes the color term 
correction in the photometric calibration is not required and, 
therefore, has not been applied. 

For each field a general catalogue was created including all 
detected sources, their magnitude in all bands and the relevant 
SExtractor parameters. Close to very bright stars the efficiency 
and reliability of source detection and photometry are signifi- 
cantly lower than the average. These regions were masked au- 
tomatically using SExtractor parameters that describe the shape 
and the size of the bright stars and sources detected in these ar- 
eas were deleted from the catalogue. Finally, visual inspection of 
the detected sources was performed on the images in order to re- 
move spurious detections and to verify the detection uniformity 
across the effective area. 

To separate stars from galaxies we used the SExtractor 
neural-network classifier fed with isophotal areas and the peak 
intensity of the source. The classifier returns a parameter named 
classstar with a value between 1 (for a point-like object) and 
(for an extended profile). Fig. Q] shows the classstar parameter 
as function of magnitude for the sources detected in one of our 
co-added image. For bright objects stars are clearly divided by 
galaxies while for fainter magnitudes the sequences spread out 
and merge. An intrinsically wide scatter in classstar is typical 
of faint galaxies. 

We defined as g alaxies all sources having classstar < 0.9 
dArnouts et alj|200lt) . In principle this choice may result in an 
overestimate of the number of galaxies at faint magnitudes but 
ensures that our sample is not biased against compact galaxies. 
Notice however that because of the relatively bright cut-off limit 
adopted for the galaxy sample selection (see Fig.Q~]and Sect. [43} 
our analysis is confined to a regime where the star/galaxy sepa- 
ration is reliable. 



4.2. Catalog crosschecks 

To check the accuracy of the astrometric and photometric cali- 
bration of our catalogs we performed a number of comparisons 
with the results of galaxy surveys and with model predictions. 

In particular one of our fields, dubbed AXAF, partially over- 
laps with the Chandra Deep Field South which has bee ob- 
served also with WFI at the ESQ2.2m teles cope during the ESO 
Imagin g Survey (EIS. lArnouts et al.ll200ll) and the COMBO-17 
survey (Wo lfet al. 112003b . In the overlapping area we identified 
~ 1200 galaxies in common with the COMBO-17 R band source 
list. The comparison of the astrometric calibration between our 
catalog and COMBO-17 ( Fig. |2]i shows a small scatter (r.m.s. 
~ 0.15 arcsec) both in right ascension and in declination. Also 
the comparison of the R band photometric calibration (Fig. |3]l 
shows a small zero point offset (~ O.Olmag) and a small disper- 
sion (rms ~ 0.03mag) for magnitudes brighter than the cut-off 
limit (fl = 21.8 mag). 

To verify the photometric calibration for all the other fields, 
we selected all bright (non saturated) stars, i.e. sources with 
classstar > 0.9 and brighter than V = 19 mag, and compared 
them to the Landolt's stars on two colour-colour (B-V, V—R,V— 
R,R - I) diagrams. The magnitudes of the Landolt's stars were 
transformed onto the WFI filter system using the correspond- 
ing color terms as reported in the instrument web pagfl We also 
checked the observed number counts as a function of colour with 
those predicted by the galactic star count model of Girard fet alj 
d2005l) TRILEGAlE for the WFI filter's set. Performing these 
tests for all fields we were able to confirm the photometric cali- 
bration to better that ~ O.lmag in all bands. 



http://www.ls.eso.org/lasilla/sciops/2p2/E2p2M/WFI/zeropoints 
http://trilegal.ster.kuleuven.be/cgi-bin/trilegal 



Botticella etal.: STRESS 



5 



I I I | I I I | i I i | r 

0.4 - 



0.3 - 




OS 



-0.2 - 



-0.4 - 

i i i I i i i I i i i I i__ 

16 18 20 22 

R c=mbo ( ma §) 

Fig. 3. Difference between the R magnitudes of galaxies in our 
source list for the AXAF field and in the COMBO- 17 catalog as 
a function of the Rcombo magnitude for magnitudes brighter than 

Rcombo — 21.8. 

4.3. Galaxy sample selection and characterization 

Defining a photometric redshift catalogue involves fulfilling two 
requirements: i) to photometrically define a sample of galaxies 
for which reliable photometric redshifts can be obtained and if) 
to characterize the redshift error distribution. 

In order to estimate the galaxy redshifts we used the spectral 
energy distribution (SED) template fitting technique. In practice, 
a grid of predicted galaxy colours as a function of redshift is 
constructed performing synthetic photometry on different galaxy 
spectral templates, either empirical or produced by spectral syn- 
thesis modeling. Then, ax 2 test of the observed vs. the predicted 
colors is used to select the best fit redshift and spectral template 
in the grid for each galaxy of the catalogue. 

The SED fitting procedure is based on the fit of the over- 
all shape of spectra and on the detection of strong spectral fea- 
tures, in particular, the 4000A, the Balmer or Lyman continuum 
breaks. Obviously the accuracy of this technique improves with 
the spectral range coverage: with B to R band photometry reli- 
able redshifts can be obtained only in the range 0.2 < z < 0.8, 
which may be extended to z ~ 1 by adding the / band photom- 
etry. Redshifts for both very nearby galaxies, which need the U 
band photometry, and more distant galaxies, which require in- 
frared (IR) filter coverage, are definitely more uncertain. 

In order to reduce the contamination from distant galaxies, 
which may be erroneously placed at low or intermediate red- 
shift, we set a cut-off limit for the galaxy catalogue at an appar- 
ent magnitude of R = 21.8 mag. This corresponds to M* B at red- 
shift z — 0.8, where M* B is a parameter of the Schechter fun ction 
which fits the galaxy luminosity function (IWolf et al. 120031) . We 
stress that for computating SN rates, we do not need the galaxy 
catalogue to be complete within a given volume. Nevertheless 
we note that at redshift z ~ 0.2 - 0.3, where the observed SN 
distribution peaks (see Sect l5.2b . the contribution to the total lu- 
minosity from galaxies fainter than R = 21.8 mag is only 25%. 



Photometric redshifts for t he selected galaxies w ere com- 
puted using the code hyper-£\ (B olzonella et al. Il2000l) running 
in the predefined redshift range < z < 0.8. As spectral tem- 
plates we tested bo th the observed spectra of Colem an, Wu & 
Weedman (CWW) (IColeman, Wu & Weedmaiill 19801) as well as 
the synthetic models based on Bruzual & Chariot (BC) library 
dBruzual & Chariot I fl993l) . Along with the redshift, the code 
provides the rest frame B-band absolute magnitude for the best 
fitting template. To further reduce the contamination from out- 
liers we set a prior on the galaxy absolute luminosity. When Mb 
is out of the range (-23, -16) mag we examine the second solu- 
tion provided by hyper-z, and if this is also inconsistent with the 
luminosity prior, we remove the galaxy from the catalog, This 
happened to 16% of the galaxies in the sample. 

Our final galaxy sample counts 43283 objects. 

The accuracy of the photometric redshifts is estimated by 
comparing them to spectroscopic redshifts, when available, and 
to the photometric redshifts obtained by the COMBO- 17 survey, 
which has a better SED sampling and used a different photo-z 
code dWolf etal.ll2003h . 

The spectroscopic redshifts for 470 galaxies of our sam- 
ple were retrieved from the Nasa/IPAC Extragalactic Database^, 
and complemented with 92 unpublished redshift measurements 
from the EDISC survey (Poggianti, private communication). We 
quantify the reliability of the photometric redshifts by measur- 
ing the fractional error, Az = (z p h - z sp )/(i + z S p) where z p h and 
z sp are the photometric and spectroscopic redshifts, respectively. 
The distribution of the fractional error, both in the case of CWW 
templates (Fig. |4]l and in that of BC templates, does not show 
systematic off-sets: for both cases a (Az) = 0.01 and cr(A) ~ 0.12 
appear consistent with our limited SED sampling. This accuracy 
is sufficient for our statistical analysis. In addition, the fraction 
of "catastrophic" outliers, galaxies with Az > 3cr, is relatively 
low, ~ 3%. 

We also tested wether there are systematic trends between 
the photometric redshift reliability and redshift range and the 
galaxy spectral type and magnitude. While there is no evidence 
of systematic trends with the magnitude (because we have se- 
lected bright galaxies), the dispersion and the fraction of outliers 
increase significantly at z < 0.2, as expected, since we lack the U 
band coverage. The dispersion of Az also increases going from 
early to late type galaxies, due to the fact that the photo-z tech- 
nique strongly relies on the strength of the Balmer break, which 
is weaker in later type galaxies. 

Similar conclusions are derived from the comparison be- 
tween photometric redshifts of galaxies in the catalogue of 
AXAF and the reds hifts for the same galaxies measured in the 
COMBO-17 survey dWolf et al. 112003b as shown in FigS 

5. The SN sample 

In this section we describe how we searched for SNe in our 
galaxy sample and obtained their redshift distribution. We first 
illustrate the detection of variable sources and the selection of 
SN candidates (Sect. 15. It then their confirmation and classifica- 
tion (Sect.lBTZii. 

5.1. SN candidate detection and selection 

In order to achieve a rapid, automatic detection of SN candidates 
we developed the S TRES S package, a collection of IRAF tasks 

7 www.web.ast.obs-mip.fr 

8 http://nedwww.ipac.caltech.edu/ 
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(zsp.-zphot.)/( 1 +zsp.) 

Fig. 4. Distribution of the fractional error Az = (z p h -z)/(l +z) in 
the comparison of our estimates z p h with spectroscopic redshift 
(thin line) and with photometric redshit of COMBO- 17 (thick 
line). 



and freely distributed software. After the standard data reduc- 
tion (see Sect. 14. Il l, the three jittered exposures of a given field 
are mapped to a common geometrical grid, properly scaled in in- 
tensity and stacked together. From this search, a suitable archive 
image (template) obtained at a different epoch is subtracted, and 
the difference image is searched for variable sources. 

After the images have been geometrically registered and 
photometrically scaled, an accurate match of the PSF of the two 
imag es is obtaine d using the mrj_phot task from the ISIS pack- 
age (jAlard 2000) which derives a convolution kernel comparing 
selected sources from the two images. The best results in the 
image subtraction are obtained when the two images have simi- 
lar seeing. For this reason we maintained an on-line archive of 
template images with different seeing. We note however that all 
the selected SN candidates were still detected even when using 
a different template image. 

Variable sources are detected as positive residuals in the dif- 
ference image using SExtractor with a source detection thresh- 
old of 3cr above the background noise and a minimum of 10 
connected pixels above the threshold. The choice of the detec- 
tion threshold is a trade-off between detecting sources as faint as 
possible and limiting the number of spurious detections. Indeed, 
the list produced by SExtractor is heavily contaminated by spu- 
rious sources due to imperfect removal of bright stars, cosmic 
rays, hot or dead pixels. In the catalogue obtained from the dif- 
ference image, false detections outnumber real variable sources 
by a factor of ~ 100. In addition, besides the SN candidates the 
detection list contains other variable sources such as fast moving 
objects, variable stars, variable AGNs. 

To clean up the detection list we used a custom-made rank- 
ing program which assigns a score to each source based on 
several parameters measured by SExtractor in the search, tem- 
plate and difference images: the class_star, the full width at half 
maximum (FWHM), the distance of the residual with respect to 



the center of the associate galaxy (if any), and apparent magni- 
tudes of the residual measured with different prescriptions. The 
scoring algorithm, tuned through extensive artificial star exper- 
iments, produces a sorted detection list. The final selection of 
trusted variable sources, about ten per image, was performed 
through visual inspection by an experienced observer. 

In particular, the residuals of fast moving objects, which 
were not completely masked by our jittering strategy, have ir- 
regular shape in the stacked image and are easily recognized by 
comparison of the individual dithered images. Conservatively, 
we classified as variable stars those sources showing a stellar 
profiles both in the search and in the template images. However, 
it is impossible to distinguish AGNs from SNe that exploded 
near galaxy nuclei by examining images taken at one epoch only. 
We labeled all variable sources closer than 0.5 arcsec to the host 
galaxy nucleus as SNAGN candidates and maintained them in 
the follow-up target list to reduce biases in the SN candidate 
selection as much as possible. The role of AGN contamination 
for those SNAGN candidates without spectroscopic observations 
will be discussed in Sect. 15.21 Finally we classified as SN can- 
didates those variables with a stellar profile in the difference im- 
age, no matter if they appeared projected on a (host) galaxy in 
the search and in the template images, or not. However, in the 
latter case we excluded the candidate from our analysis, since 
we concentrate on SNe occurring in the galaxies within our spe- 
cific sample. 

To deal effectively with the large number of epochs and can- 
didates, we develop ed a MySQI^ database with a web interface 
(Riell o et al. 1 12003) to easily access the information remotely 
(e.g. during observing runs). In particular the database stores: the 
monitored fields, pointing coordinates, finding charts, log of ob- 
servations (with seeing, exposure time, photometric zero point), 
all discovered variable sources, their identification chart and rel- 
evant parameters (sky coordinates, epoch, magnitude, stellarity 
index, etc.). 

At the end of our SN search programme we reviewed all vari- 
able sources recorded in our database to obtain a final classifi- 
cation based on all information gathered during our observing 
programme. A search engine allows us to identify multiple de- 
tections of the same source in different epochs and filters. This 
was used for the identification of AGNs which in general show 
long term, irregular variability. Indeed, the relatively long dura- 
tion of our search makes this approach fairly efficient for AGN 
removal. 

5.2. SN candidate confirmation and classification 

Spectroscopic observations were planned for all SN and SNAGN 
candidates but could not be completed due to limited telescope 
time and, in some cases, poor weather conditions. In general a 
higher priority was given to candidates flagged as SN. Spectrum 
extraction and calibration were performed using standard IRAF 
packages. Special care was devoted to removal of background 
contamination through a detailed analysis of the intensity profile 
along the slit. For the nuclear candidates it is impossible to sepa- 
rate the variable source from the background using this informa- 
tion only. In such cases, first the sky emissions was subtracted, 
then the spectrum was extracted by integrating all the light along 
the slit. A best fit galaxy spectral template extracted from the 
STSci Database of UV-Optical Spectra of Nearby Quiescent and 



9 MySQL is an open source database released under the GNU 
General Public License (GPL), http : / /www . mysql . com/ 
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Table 3. Spectroscopically confirmed SNe 



Design. R.A.(J2000.0) Dec.(J2000.0) type z field band reference 



1999ey ff 


00:58:03.43 


-27:40:31.1 


Iln 


0.093 


J1888 


V 


IAUC7310 


1999gt 


03:32:11.57 


-28:06:16.2 


la 


0.274 


AXAF 


V 


IAUC 7346 


1999gu 


03:33:00.20 


-27:51:42.7 


11 


0.147 


AXAF 


V 


IAUC 7346 


2000fc 


00:58:33.55 


-27:46:40.1 


la 


0.420 


J1888 


V 


IAUC 7537 


2000fp 


05:38:01.27 


-23:46:34.1 


11 


0.300 


EisC 


V 


IAUC 7549 


2001bd 


19:13:10.94 


-64:17:07.8 


II 


0.096 


Field2 


V 


IAUC 7615 


200 lgf 


22:04:21.27 


-18:21:56.8 


la 


0.200 


22Z3 


V 


IAUC 7762 


2001gg 


00:45:20.83 


-29:45:12.2 


II 


0.610 


EisB 


V 


IAUC 7762 


200 lgh 


00:57:03.63 


-27:42:32.9 


II 


0.160 


J1888 


V 


IAUC 7762 


200 lgi* 


04:28:07.06 


-36:21:45.2 


la 


0.200 


EF0427 


V 


IAUC 7762 


2001gj 


05:38:07.90 


-23:42:34.4 


II 


0.270 


EisC 


V 


IAUC 7762 


200 lio 


00:45:53.64 


-29:26:11.7 


la 


0.190 


EisB 


V 


IAUC 7780 


200 lip 


03:31:13.03 


-27:50:55.5 


la 


0.540 


AXAF 


V 


IAUC 7780 


2002cl 


13:44:09.94 


-00:12:57.8 


Ic 


0.072 


13Z3 


V 


IAUC 7885 


2002cm 


13:52:03.67 


-11:43:08.5 


11 


0.087 


WH23 


V 


IAUC 7885 


2002cn 


13:53:21.24 


-12:15:29.5 


la 


0.302 


WHhz2 


V 


IAUC 7885 


2002co* 


14:10.53.04 


-11:45:25.0 


11 


0.318 


WH27 


V 


IAUC 7885 


2002du 


13:53:18.28 


-11:37:28.8 


II 


0.210 


WH23 


V 


IAUC 7929 


2004ae 


04:28:17.89 


-36:18:55.0 


II 


0.480 


EF0427 


R 


IAUC 8296 


2004af 


05:38:03.91 


-23:59:00.2 


Ic 


0.056 


EisC 


R 


IAUC 8296 


2004ag 


09:51:01.49 


-20:50:37.5 


11 


0.362 


EisD 


R 


IAUC 8296 


2004ah 


10:45:47.41 


-00:06:58.1 


la 


0.480 


10Z2 


R 


IAUC 8296 


2004ai 


13:54:26.09 


-12:41:15.9 


Ic 


0.590 


WHhz2 


R 


IAUC 8296 


2004aj* 


14:20:37.58 


-12:24:14.4 


la 


0.247 


WH31 


R 


IAUC 8296 


2004cd* 


13:44:50.70 


-00:03:48.3 


la 


0.500 


13Z3 


R 


IAUC 8352 


2004ce 


13:54:45.93 


-12:14:36.3 


la 


0.465 


WHhz2 


R 


IAUC 8352 


2004cf 


14:11:05.77 


-11:44:09.4 


Ib/c 


0.247 


WH27 


R 


IAUC 8352 


2004cg 


14:11:23.71 


-12:01:08.4 


II 


0.264 


WH27 


R 


IAUC 8352 


2004dl 


14:11:11.50 


-12:02:36.9 


la 


0.250 


WH27 


R 


IAUC 8377 


2004dm 


22:43:21.36 


-40:19:46.7 


lb 


0.225 


EisAl 


R 


IAUC 8377 


2005cq* 


09:52:00.48 


-20:43:26.5 


la 


0.310 


EisD 


R 


IAUC 8551 



# - Not included in the rate computation because the host galaxy has R > 21.8. 



Active Galaxies was then subtracted from the candidate spec- 
trum. 

The classification of SN candidate spectra was facilitated 
by the use of Passpartou, a software package that automati- 
cally classifies the candidate spectrum after comparison with 
all the SN spectra in th e archive of the Padova SN group 
dHarutvunvan et alj|2005l) . 

We could obtain direct spectroscopic observations only for 
38 candidates (~ 40% of the total): 31 were confirmed to be 
SNe whereas 7 turned out to be AGN. In this paper we use only 
the information on SN type and redshift of the host galaxy while 
in a future paper we will give more details on the spectroscopic 
data analysis. 

The properties of the 3 1 confirmed SNe are summarized in 
Table [3] 25 have been observed in our galaxy sample (9 SNe la 
and 16 SNe CC) whereas 6 SNe discovered in galaxies with R > 
21.8 mag were excluded from the following analysis. The SN 
redshifts range from z = 0.056 to z = 0.61. 

As mentioned in Sect[3] we carried out a complementary ob- 
serving programme to estimate redshift and spectral type of SN 
host galaxies and to better constrain contamination of AGNs in 
the SNAGN candidate sample. We remark that the higher prior- 
ity given to SN candidates for immediate follow-up introduces 
a bias against SNAGN candidates. Such bias is minimized by 
this late follow-up which mostly targets the host galaxies of the 
SNAGN candidates. We found that, out of 42, 22 host galaxies 
showed an active nucleus. All the candidates occurring in these 
galaxies were removed from the candidate list. Actually we have 

10 http : //www . stsci . edu/f tp/catalogs/nearby_gal/ sed . html 



to notice that with our approach a SN occuring in a galaxy with 
an active nucleus will be discharged if only late spectroscopy is 
available. In Table|4]we report coordinates, redshifts, search fil- 
ter and classifications of the remaining 20 SN candidates with 
host galaxy spectroscopy. 

Even after this complementary analysis, 13 SN and 28 
SNAGN candidates still remained without direct or host galaxy 
spectroscopy. The list of these remaining candidates is reported 
in Tab. 

For all these candidates the redshift of the host galaxy was 
estimated using the photo-z technique. To avoid a possible con- 
tamination by the SN candidate in measuring the host galaxy 
colours we produced deep images stacking only those frames 
obtained 6-12 months earlier than the epoch of discovery of the 
candidate. 

In summary the analysis presented here is based on 86 ob- 
jects: 25 spectroscopically confirmed SNe, 33 SN candidates and 
28 SNAGN candidates. The redshift distribution of candidates 
peaks at z ~ 0.2 - 0.3 and extends up to z ~ 0.6 as shown in Fig 
[5] Combining information from the long term variability history, 
direct and host galaxy spectroscopy we found that ~ 50% of 
the variable sources originally classified SNAGN were actually 
AGN. Hereafter, in the analysis of SN statistics, we assigned a 
weight of 0.5 to all SNAGN candidates with no spectroscopic 
observations. The unclassified SN and SNAGN candidates were 
distributed among type la and CC based on the observed frac- 
tions in each redshift bin for the SN subsample with spectro- 
scopic confirmation. It turns out that, with respect to the total 
number, SNe la are 27 ± 18%, 47 + 21% and 63 + 36% at redshift 
z = 0.1, 0.3, 0.5 respectively. We quantify the effect of misclassi- 
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Table 4. SN candidates with host galaxy spectra Table 5. Candidates without direct or host galaxy spectroscopy 



Design. 


R.A.(J2000.0) 


Dec.(J2000.0) 




z 


band 


03Z3.B 


03 


39 


27 


76 


-00:01:08.8 





240 


V 


03Z3.F 


03 


40 


36 


80 


-00:06:17.8 





247 


V 


03Z3.H 


03 


40 


09 


43 


-00:10:20.6 





282 


V 


03Z3J 


03 


40 


37 


11 


-00:05:22.5 





180 


V 


10Z2_B 


10 


45 


42 


76 


00:00:27.9 





295 


V 


10Z2JD 


10 


47 


40 


94 


-00:13:52.6 





350 


V 


EF0427.G 


04 


30 


21 


93 


-36:06:49.9 





359 


V 


EF0427.I 


04 


30 


10 


73 


-36:25:47.4 





420 


V 


EF0427J 


04 


29 


20 


47 


-36:32:16.7 





185 


V 


EF0427.O 


04 


28 


10 


56 


-36:13:54.4 





450 


V 


EF0427.Q 


04 


28 


17 


09 


-36:09:54.9 





394 


V 


EisAl_B 


22 


44 


01 


94 


-40:03:21.7 





219 


V 


EisAl.C 


22 


44 


25 


26 


-40:20:01.9 





215 


V 


EisAU 


22 


44 


02 


82 


-40:02:05.2 





180 


V 


EisAl.Q 


22 


42 


37 


87 


-39:57:54.4 





380 


V 


EisAl_R 


22 


42 


13 


27 


-40:03:19.9 





530 


V 


EisC_A 


05 


38 


45 


67 


-23:44:41.0 





643 


V 


EisC.C 


05 


38 


06 


55 


-23:36:57.1 





133 


V 


EisC.G 


05 


37 


52 


11 


-23:38:14.4 





156 


V 


WH27_B 


14 


11 


02 


77 


-11:49:06.2 





485 


V 


WHhz2.A 


13 


53 


58 


73 


-12:43:04.8 





444 


V 


WHhz2.B 


13 


54 


25 


27 


-12:35:42.9 





322 


V 




red shift. 

Fig. 5. The redshift distribution of SN candidates with (la and 
CC) and without spectroscopic confirmation (SN and SNAGN). 



fixation for the candidates without spectroscopic observations by 
performing Monte Carlo simulations as illustrated in Sect. 17.21 

6. The control time of the search 

The time interval during which a SN occurring in a given galaxy 
can be detected (the control time) depends mainly on the SN 
detection efficiency (Sect. 16.11 ). the SN light curve (Sect. I6.21 > 
and the amount of dust extinction (Sect. l6~4b . In this section we 
describe the different ingredients needed to calculate the control 
time of each single observation and the recipe to obtain the total 
control time of our galaxy sample (Sect. l6~5l ). 

6.1. The search detection efficiency 

The SN detection efficiency depends mainly on the survey strat- 
egy, instrumental set up and observing conditions. For a given 



Design. 


R.A.(J2000.0) 


Dec.(J2000.0) 


z 


band 


class. 


03Z3.Q 


03:38:26.39 


-00:04:53.1 


0.30 


V 


snagn 


10Z2.A 


00:47:06.51 


00:00:39.6 


0.52 


R 


snagn 


13Z3.R 


13:44:01.94 


00:00:03.9 


0.35 


R 


snagn 


13Z3.X 


13:45:22.55 


00:03:11.2 


0.18 


R 


snagn 


13Z3.AB 


13:44:31.26 


00:00:57.5 


41 


R 


sn 


13Z3.AC 


13:45:11.26 


-00:12:27.4 


0.43 


R 


snagn 


22Z3.I 


22:06:03.26 


-18:28:10.8 


0.42 


R 


sn 


22Z3J 


22:04:42.12 


-18:36:38.3 


0.19 


R 


sn 


22Z3.K 


22:04:35.87 


-18:45:45.0 


0.42 


R 


snagn 


22Z3.N 


22:04:49.66 


-18:45:27.3 


0.14 


R 


snagn 


22Z3.R 


22:04:03.87 


-18:35:55.2 


0.22 


R 


sn 


AXAF_B 


03:32:31.16 


-28:04:43.8 


0.13 


V 


sn 


AXAFJD 


03:31:28.98 


-28:10:26.1 


0.42 


V 


snagn 


AXAF_F 


03:33:05.29 


-27:54:09.2 


0.17 


V 


snagn 


AXAF.G 


03:32:39.23 


-27:42:57.5 


0.47 


R 


sn 


AXAFJSf 


03:32:47.68 


-28:07:41.1 


0.18 


R 


sn 


EF0427_B 


04:29:07.28 


-36:28:01.8 


0.50 


V 


sn 


EF0427_F 


04:29:49.31 


-36:33:55.2 


0.24 


V 


snagn 


EF0427JV1 


04:27:57.15 


-36:30:52.9 


0.74 


V 


snagn 


EisAl_E 


22:43:44.86 


-40:05:50.6 


0.17 


V 


snagn 


EisAU 


22:44:48.37 


-40:03:42.3 


0.17 


V 


snagn 


EisAl.O 


22:42:01.73 


-40:21:17.8 


0.25 


V 


snagn 


EisAl_AC 


22:44:21.82 


-40:08:21.4 


0.52 


R 


sn 


EisB_B 


00:45:09.35 


-29:28:33.2 


0.13 


V 


snagn 


EisB_P 


00:44:14.57 


-29:39:08.1 


0.17 


V 


snagn 


EisB_S 


00:45:51.29 


-29:52:10.4 


0.43 


V 


snagn 


EisC_H 


05:39:28.83 


-23:38:05.7 


0.36 


V 


sn 


EisCJ 


05:38:01.65 


-23:51:40.7 


0.20 


V 


snagn 


EisCJC 


05:38:28.78 


-23:46:29.4 


0.80 


V 


snagn 


EisCLO 


05:39:28.51 


-24:00:05.4 


0.80 


V 


snagn 


EisC_Q 


05:37:50.49 


-23:36:55.2 


0.14 


R 


sn 


EisD_F 


09:51:12.26 


-20:56:14.0 


0.16 


R 


snagn 


Field2.C 


19:11:16.62 


-64:17:32.3 


0.16 


V 


snagn 


Field2_E 


19:10:35.43 


-64:21:36.0 


0.28 


R 


sn 


J1888J 


00:58:41.35 


-27:50:38.1 


0.30 


V 


sn 


J1888.M 


00:57:05.34 


-27:45:57.7 


0.63 


V 


snagn 


WH23.A 


13:53:18.79 


-11:40:29.7 


0.16 


V 


snagn 


WH23.E 


13:52:41.23 


-11:51:24.7 


0.21 


R 


snagn 


WH27.N 


14:11:47.81 


-11:35:25.6 


0.44 


R 


snagn 


WHhz2_Q 


13:54:31.46 


-12:18:22.5 


0.67 


R 


snagn 


WHhz2_T 


13:55:08.90 


-12:18:38.3 


0.23 


R 


snagn 



instrumental set up and exposure time, the detection efficiency 
curve as a function of the SN candidate apparent magnitude was 
estimated via Monte Carlo simulations. 

In each simulation artificial point-like sources of different 
magnitudes were added to an image which was then searched for 
variable sources using the same software as in the actual search. 
The detection efficiency at a given magnitude is computed as 
the ratio between the number of recovered and injected artificial 
sources. 

Artificial point-like sources were generated using the IRAF 
Daophot package with the PSF determined from a number of 
isolated stars on each chip of each image. The artificial stars 
were then placed on the galaxies (one source per galaxy) in a 
random position with respect to the center, but following the the 
luminosity profile. Detailed simulations for a few selected fields 
allowed us to probe the most relevant parameters affecting the 
detection efficiency: the search image characteristics (observing 
conditions and effective search area), the search process (choice 
of template image and selection criteria) and the SN position in 
the host galaxy. 
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23.4 
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Table 6. Adopted B band, absolute magnitude at maximum for 
different SN types. 



22.8 - 



22.6 



seeing 



Fig. 6. The magnitude at which the detection efficiency is 50%, 
m$o%, in R band as function of seeing. 



In particular we found that, for a given seeing, differences 
in the sky transparency do not affect the shape of the efficiency 
curve. The magnitude at which the efficiency is 50% (m.50%) be- 
comes fainter when the sky transparency improves, and a vari- 
ation of photometric zero point translates directly into a change 
in m 50 %. 

After accounting for the differences in sky transparency, a 
strong dependence of the detection efficiency on seeing is still 
present (Pig. [6j>. We found that mso% is ~ 1 mag brighter when 
the seeing degrades from 0.7 arcsec to 1.2 arcsec, as shown in 
Fig-El If the seeing is worse than 1.3-1.4 arcsec then 11150% de- 
creases more rapidly. 

The detection efficiency for different fields is very similar ex- 
cept for minor variations due to different effective search areas 
caused by the presence different numbers of bright stars (which 
were masked). This is also the reason why the detection effi- 
ciency is always lower than 100%, even for bright sources. 

In order to reduce the overhead of determining the detec- 
tion efficiency for each individual observations of a given field 
we performed an extensive set of Monte Carlo simulations for 
only one field, covering a wide range of source magnitudes and 
seeing, thereby generating a set of template efficiency curves. 
Then, for each individual observation of the others fields, we 
performed only few simulations to determine the shift to apply 
to the appropriate template curve. 

Possible spurious effects caused by the choice of the tem- 
plate image were tested by performing the same analysis based 
on three other images with different seeing. We found that, in the 
range 21 -23 mag, 95% of the sources are recovered in all cases, 
while, at fainter magnitudes, this fraction quickly decreases. 

As we mentioned in Sect. 15. II the final selection of SN candi- 
dates was based on visual inspection. We verified that the detec- 
tion efficiency is independent of the observer that carried out the 
visual inspection by comparing the list of candidates produced 
by different observers in several simulations. Finally, we found 
that the detection efficiency depends very little on the position of 



SN type M B . 



<A B > 



Ref. 



Ia 
Ib/c 
IIP 
I1L 



-19.37 
-18.07 
-16.98 
-18.17 



0.47 
0.95 
1.00 
0.53 



0.4 
0.7 
0.7* 
0.7* 



Wang et al. (2006) 
Richardson et al. (2006) 
Richardson et al. (2002) 



#- Nugent et al. (2006) 

the source with respect to the galaxy and with galaxy brightness 
except for nuclear candidates in very bright galaxies. 

This analysis of the detection efficiency was carried out both 
on R- and V-band images with very similar results after consid- 
ering the zero-point differences. 

6.2. SN light curves 

The observed light curve in the band F of a given SN type at 
redshift z is given by: 



mf(t,z) = M B N (0) + AM B N (t') + Kf F (z,t') 



+Af+A B+f iCz) 



(1) 



where: 



- M^ N (0) is the SN absolute magnitude at maximum in the B 
band, 

- AM™(f') describes the light curve relative to the maximum, 

- K B ^(z, t') is the K-correction from the B to the F band, (V or 

R), 

- A^. and A h B are the Galactic and host galaxy absorption re- 
spectively, 

- fi(z) is the distance modulus for the redshift z in the adopted 
cosmology. 

The light curve is translated into the observer frame by ap- 
plying the time dilution effect to the rest frame SN phase t' = 
f/(l + z). Four basic SN types were considered: Ia, IIL, IIP, Ib/c. 
W e adopted the same tem plates for the SN B-band light curves as 
in lCappellaro et alj d 19971) and account for the observed disper- 
sion of SN absolute magnitudes at maximum light by assuming 
a normal distribution with the mean Mb ,0 and cr listed in Tab. [6] 
In this table we also report the average absorption < Ag > mea- 
sured in a sample of nearby SNe. 

Finally, Ga l actic extinction is taken into account using the 
ISchlegel et al. I d 1998b extinction maps, whereas the correction 
of the host galaxy extinction is estimated through a statistical 
approach, described in Sect. [ 



6.3. K correction 

For a given filter, ^-correction is the difference between the 
magnitude measured in the observer and in the source rest frame 
which depends on the redshift and SED of the source. In princi- 
ple, if the source redshift is known it is possible to design a filter 
properly tuned to minimise the need for a ^-correction. In most 
practical cases, however, sources are distributed in a wide range 
of redshifts and this approach is inapplicable. 

A detailed discussion of the formalism and application of 
.g-correction to the case of type Ia SNe can be found in 
iNugent etal.l (120021) . SN Ia have received a special care as their 
use as cosmological tools requires high precision photometry. 
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Actually, for computing the control time, we can tolerate much 
larger uncertainties. However, we need to obtain the K correc- 
tion also for type Ib/c and type II SNe, which have not been 
published so far. 

In general, there are three main problems for an accurate esti- 
mate of the /f-correction of SNe: i) the SED rapidly evolves with 
time it) SNe, in particular type Ib/c and type II, show large SED 
diversities, Hi) the UV spectral coverage, which is needed to 
compute optical Zf-correction at high redshift, is available only 
for few SNe, and typically only near maximum light. 

We addressed these issues trying to use all available infor- 
mations on SN SED, in particular: 

1. we selected a sample of spectral sequences of different SN 
types with good S/N retrieved from the Asiago-Padova SN 
archive. The sample includes the type la SNe 1990N, 1991T, 
1991bg, 1992A, 1994D, 1999ee, 2002bo, the type II SNe 
1979C, 1987 A, 1992H, 1995G, 1995 AD, 1996W, 1999el, 
2002GD, 2003G and the type Ib/c SNe 19901, 1997B, 
1997X, 1998bw. 

2. UV spectral coverage was secured retrieving SN spec- 
tra from the ULDA SN spectra archive (SN la 1990N 
phase -9, 1981B phase 0) of th e International Ultraviolet 
Explorer ( Cappellar o et alJl 1995b and from the HST archive 
(SN la 1992A at phase +5, +45, SN Ib/c 19941 phase +7d, 
SN 1999em phase 5d) . These were combined with optical 
and IR spectra to obtain the SED in a wide range of wave- 
lengths. 

3. we include d models for sta ndard SNIa at phase -7, 0, +15 d 
provided bv fMazzali I d2000l) and mod els of the SNI I 1999em 
at phase 0,+9,+l l,+25 d, courtesy of iBaron et all d2004l) . 

The observed and model spectra were shifted at different red- 
shifts and synthetic photometry in the required bands was ob- 
tained using the IRAF package synphot. Eventually, we derived 
the Kf^ and Kf£ corrections for each SN type as a function of 
light curve phase with a step of 0.05 in redshift. Examples of the 
measured K correction as a function of the SN phase are shown 
in on line Figs.[TlfT8l 

Uncertainties due to variances of the SN spectra, uncertain 
extinction correction and/or incomplete temporal/spectral cover- 
age can be minimized by choosing a proper filter combination. 
Indeed as it can be seen in on-line Fig.[14l for the average red- 
shift of our search, z ~ 0.3, K B y is almost independent on the 
SN type and phase. 

6.4. Host galaxy extinction 

The host galaxy extinction correction is the most uncertain in- 
gredient in the estimate of the SN rate. In the local Universe 
Capp ellaro et al.l 0999) relied on an empirical correction as a 
function of the host galaxy morphological type and inclination. 
This approach cannot be applied to galaxies of our sample be- 
cause the relevant data are not available. Hence we resorted to a 
statistical approach, based on the modelling of SN and extinction 

distribution in galaxies. 

In short, following the method described in Riello & Patat I 
(2005), we performed a number of Monte Carlo simulations 
where artificial SNe were generated with a predefined spatial 
distribution function and were seen from uniformly distributed 
lines of sight. Integrating the dust column density along the line 
of sight for each SN we derived the total optical depth and the 
appropriate extinction was applied to the SN template spectrum. 
Repeating a number of simulations we obtained the expected 
distribution of SN absorptions. 



Table 7. Parameters of our modelling of the host extinction for 
the four SN types considered in the simulati ons. The last column 
of the table provides the reference from iRiello & Patat I d2005l) 
where the meaning of each parameter is explained in detail. 



Parameter 


la 


Ib/c 


IIP 


Dust 


Ref. 


B/T 


0.5 










§2.1 




1.8 








Eq. (1) 


n 


8 








§2.1.1 




4.0 


4.0 


4.0 


4.0 


Eq. (2) 




0.35 


0.35 


0.35 


0.14 


Eq. (2) 


n 


6 


6 


6 


6 


§ 2.1.2 


m 


6 


6 


6 


6 


§2.1.2 


N u 





2 


2 


2 


Eq. (3) 


w 




0.2 


0.2 


0.4 


Eq. (3) 


P 




20 


20 


20 


Eq. (3) 



With respect to IRiello & Patat I d2005l) . we included in the 
modelling also CC SNe. Such extension was straightforward: we 
added spectral templates for CC SNe, and adopted a reasonable 
spatial distribution for the simulated CC SNe. Since CC SNe 
are thought to originate from massive progenitors, the spatial 
distribution is assumed to be concentrated in dust-rich regions 
along the spiral arms. 

The parameters for SN and dust spatial distributions used 
in the simulations are summarized in Tab. [7] The top rows (1- 
3) of the Table show the value of the parameters that describe 
the bulge component, that is the bulge to total luminosity ratio 
{B/T), the scale and the truncation radius of SN distribution (r h 
and n respectively). The bottom half of the table (rows 4-10) 
lists the parameters for the disk component of both SN and dust, 
namely the scale-length and scale-height of the disc {r cl ,z d ), the 
respective truncation radii (n,m) and t he spiral arm perturb ation 
parameters. The reader should refer to Riel lo & Patat I d2005l) for 
a more detailed description of N a ,w,p parameters. 

The two other parameters governing the dust properties are 
the total to selective extinction ratio Ry and the total optical 
depth along the galaxy rotation axis r(0), which provides a 
convenient parameterization for the total amount of dust in the 
galaxy. For the former we adopted the canonical value Ry = 3.1. 
For the latter we considered two scenarios: a standard extinction 
scenario with r(0) =1.0 for both SN types and a high extinction 
scenario with t(0) = 5.0 only for CC SNe. 

For each SN type we ran a set of simulations with 10 5 artifi- 
cial SNe covering the redshift range z = 0.05 - 0.80 with a step 
in redshift of 0.05. In order to speed-up the computation, rather 
th an considering a libra ry of SN spectra at different phases as 
in IRiello & Patat I ([2005), we used only a single spectrum close 
to maximum light for each SN type. This is a reasonable ap- 
proximation since the evolution of the SN SED with time has 
a negligible effect on the absorption (at least for the SN phases 
that are relevant for our search). 

Once a simulation for a given parameter set was completed, 
we derived the observed distribution of SN absorptions, and used 
it to compute the control time (cf. Sect. 16. 5b . For t(0) = 1.0, the 
average absorption from the modeling is < A h B >= 0.5 mag for 
SNe la and < A h B >= 0.7 mag for CC SNe, after excluding SNe 
with absorption larger than A h B — 3 mag (about 8% of SNe la 
and 12% of CC SNe in this model, which are severely biased 
against in typical optical searches). We note that these average 
absorptions are consistent with the observed ones, reported in 
Tab. [6] For t(0) = 5.0 the average model absorption is < A h B >= 
0.8 mag for CC SNe with A h B < 3 mag (about 50% of the CC 
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SNe) which is still consistent with the observed value in the local 
Universe. 

Our correction for dust extinctions is meant to provide a rea- 
sonable estimate of this effect. A more realistic model should 
consider the different amount and distribution of the dust in each 
galaxy type. This is far beyond the scope of the current analysis. 



(0.056-0.61), we can obtain an observational constraint on evo- 
lution of the rate by analysing the redshift distribution of the 

events. 

In analogy with iPain et aTl d2002h. and following analysis 
dDahlen et al. ll2004tlGal-Yam & Maoz 1120041) . we adopt the fol- 
lowing power law for the reshift dependence of the rate: 



6.5. Control time calculation 

We now detail how all the ingredients described in the previous 
sections have been combined to compute the control time of our 
galaxy sample. First, for a given SN type and filter F, the control 
time of the /-observation of the y'-galaxy was computed as: 



CT 



SN,F 



J" iv N ' F (m) ej(m) dm 



(2) 



where r^ N F (m) is the time spent by the SN in the magnitude 

range m and m + dm and ef(m) is the detection efficiency. More 
specifically, we convolved the distribution of the absolute mag- 
nitude at maximum (M B N (0)) and of the absorption due to host 
galaxy extinction (A B ), so as to determine the distribution of the 
quantity M| N (0) + A B appearing in Eq.Q] and computed iv N ' F (m) 
as a weighted average of the individual times over this combined 
distribution. 

Then, the total control time Cr| N F of the y'-galaxy was com- 
puted by summing the contribution of individual observations. 
If the temporal interval elapsed since the previous observation is 
longer than the control time, that contribution is equal to the con- 
trol time of the observation, otherwise it is equal to the interv al 
of time between the two observations (ICappellaro et al.|[l999l) . 

The total control time of the galaxy sample is obtained as the 

SN,F 

B band luminosity weighted (CT • ) average of the individual 
galaxies. 

Since we merged all CC subtypes, including type Ib/c IIP 
and IIL, the control time for CC SNe is computed as follows: 



-Ib/c 



CTj =fm/cCT j ' +fuLCTj + f llP CTj 



(3) 



where the relative fractions of the different CC subtypes 
is assumed to be constant with redshift and equal to that ob- 
served in the local Univ erse, namely 20% of Ib/c and 80% of II 
( ICappellaro et alJI 1999b. out of which 35% are IIL and 65% are 
IIP events dRichardson et alj|2002l) . 

In order to illustrate the role of the host galaxy extinction, we 
calculated the CC SN rate by adopting the control time both for 
a standard (t(0) = 1.0), and for a high (t(0) = 5.0) extinction 
scenario. 



7. SN rates 

In this section we describe our approach to estimate SN rate, 
present our results (Sect. 17.11 ) and discuss their uncertainties 
(Sect. 17.2b . We also compare our measurements to those ob- 
tained in the local Universe (Sect. 1731 1. 



7.1. SN rate estimate 

The SN rate at redshift z is given by the ratio between the number 
of observed SNe and the control time of the monitored galaxies 
at that redshift. Since our SN sample spans a wide redshift range 



,SN,F/ 



1+Z 

1 + z 



(4) 



where r(z) is the rate at the mean redshift of the search, z, and 
a is the evolution index. Writing the rate evolution in terms of its 
value at the average redshift of the search reduces the correlation 
between the two free parameters in Eq. ©■ 

The reference redshift z of the search is computed as the 
weighted average of the galaxy redshifts with weights given by 
the respective control times: 



•SN,F 



-SN,F 

Z 



r M zjC Tj 



(5) 



We obtained: 



? Ia = 30 +014 

-f-C _ n ? i +0.08 
i. -u.z.i_ (U)g 



Clearly, the lower z for CC SNe is due to their being, on av- 
erage, intrinsically fainter than SNe la; hence their control time 
at a given redshift is shorter. 

For a given SN type and filter F, the number of expected SNe 
is given by the expression: 



N' 



SN,F/ 



-SN,F 



(1 + z) 



(6) 



where the sum is extended over the n galaxies at redshift z 
and the factor 1+z corrects the rate to the rest frame. 

We compared the observed SN distribution, binned in red- 
shift, to the expected distribution performing a Maximum 
Likelihood Estimate to derive the best fit values of the param- 
eters r(z) and a, independently for SNe la and CC SNe. 

Our candidates come from V and R search programs, with 
no overlap. The rates at the reference redshift based on V and R 
candidates differ by less than 20% both for type la and CC SNe, 
which is well below the statistical error. Therefore, V and R SN 
candidates were combined togheter to improve the statistics. 

For a standard extinction, the best fit values and statistical 
uncertainties for r(z) and a are: 



- r Ia (0.30) = 0.22+° ^ h% n SNu, 



„cc 



(0.21) 



0.08 "70 ' 

0.82!^ /^SNu, 



= 4.4 



+3.6 
-4.0 



„CC 



■ 7 5 +2 - 8 



Fig. [T9l shows the 1,2,3 cr confidence levels for the two pa- 
rameters describing the evolution of SN rate, whereas Fig. [20] 
shows the best fit of the observed redshift distribution separately 
for SNe la, CC SNe and then all together. 

The assumption of a particular host galaxy extinction sce- 
nario is one of the most important sources of uncertainty for SN 
rate measurements. We tested the dependence of the results on 
this effect by performing the fit under extreme assumptions: no 
extinction and a high extinction scenario, the latter only for CC 
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Fig. 19. Confidence levels for the parameters describing the red- 
shift evolution of the SN rates, from the maximum-likelihood 
test. 

SNe.The results are reported in Tab. [8] With respect to the stan- 
dard case, the type la and CC SN rates decrese respectively by a 
factor of 1.7 and 2 when no extinction is adopted. The effect is 
smaller for the SN la rate because SNe la occur, on average, in 
environments with a smaller amount of dust. If a high extinction 
correction is adopted, the CC rate increases by a factor of 2: this 
can be regarded as a solid upper limit. 

7.2. Systematic Uncertainties 

When the size of the SN sample is large enough, the statistical 
uncertainty is relatively small, and it becomes important to ob- 
tain an accurate estimate of systematic uncertainties. These were 
evaluated via Monte Carlo simulations. We identified the follow- 
ing sources of systematic uncertainties, for which we specify the 
distribution adopted for the simulations: 

- re-distribution of unclassified SNe in SN types (Gaussian 
distribution with cr from the statistical error, as discussed in 
Sect. ED), 



Fig. 20. Observed (points with statistical error bars) and ex- 
pected (solid line) redshift distribution of SN events. 

Table 8. The parameters of the rate evolution estimated for dif- 
ferent extinction scenarios. 



SN type z Extinction r(z) [/j, SNu] a 



la 


0.30 


none 


13+ 006 


7 7+3.6 
^-3.8 






standard 


22+ 010 


A 4+3.6 


CC 


0.21 


none 


0.42+°;* 


5 4+ 2 ' 8 

J -+-3.4 






standard 


0.82+J^ 


7.5+1* 






high 




7.8+1* 



- AGN contamination factor (Gaussian distribution with an av- 
erage value of 0.5 and cr = 0.25), 

- photometric redshifts (the number counts in each bin of the 
galaxy redshift distribution are left to vary by +50% with an 
uniform distribution), 

- detection efficiency (Gaussian distribution with cr = 0.1 
mag), 

- SNe absolute magnitude (Gaussian distribution with cr = 0.1 

mag), 

- standard extinction correction (uniform distribution with 
+50% variance with respect to the standard case). 

The fit described in the previous section was performed for 
each simulation, to yield a distribution of solutions for the r(z) 
and a. The results are shown in Tab. [9] for each source of sys- 
tematic uncertainty, as derived by varying only one parameter at 
the time. The total systematic error has instead been derived by 
varying all parameters simultaneously. 

Not surprinsingly, the major uncertainty is due to the lack 
of the spectroscopic classification for a large fraction of the SN 
candidates. For the CC SN rate also the estimate of the detec- 
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Table 9. Summary of systematic uncertainties 



SNIa 



SNCC 



Source item 


cr r 


cr a 


cr r 


<T„ 


SN types distribution 


+0.08 


+3.4 


+0.18 


+3.4 




-0.08 


-3.4 


-0.20 


-3.5 


AGN fraction 


+0.02 


+0.2 


+0.05 


+0.2 




-0.01 


-0.2 


-0.05 


-0.1 


Photometric redshifts 


+0.04 


+0.7 


+0.13 


+0.6 




-0.03 


-0.7 


-0.11 


-0.8 


Detection efficiency 


+0.02 


+0.6 


+0.21 


+0.6 




-0.03 


-0.6 


-0.11 


-0.8 


Standard extinction 


+0.07 


+ 1.9 


+0.26 


+0.7 




-0.05 


-1.4 


-0.22 


-0.9 


SN absolute mag 


+0.02 


+0.7 


+0.06 


+0.4 




-0.02 


-0.6 


-0.06 


-0.4 


All systematic errors 


+0.16 


+3.3 


+0.30 


+3.6 




-0.14 


-3.5 


-0.26 


-3.7 


Statistical errors 


+0.10 


+3.6 


+0.31 


+2.8 




-0.08 


-4.0 


-0.24 


-3.2 



tion efficiency and the dust extinction correction are important 
sources of uncertainty. Even for our relatively small SN sam- 
ple, the statistical and systematic uncertainties are comparable. 
Since, due to a growing number of detected SNe, the statistical 
uncertainty will decrease in the future, the systematic errors will 
soon dominate the overall uncertainty. A special care should then 
be devoted to reduce the systematic effects, in particular secur- 
ing spectroscopic follow up of all the candidates, and obtaining 
more precise information on the extinction of the galaxy sample. 

7.3. Comparison with the local SN rates 

In order to measure SN rate in a galaxy sample one needs to 
specify how the rate scales with a physical parameter propor- 
tional to t he stellar c ontent of each galaxy. In the early '70 
iTammanl d 19701: 11974 and afterwards ICappellaro et "ail d 1993b 
showed that the SN rates scale with the galaxy B band luminosity 
and from then on the SN rates, in the local Universe, have been 
measured in SNu. For a direct comparison of the rate at interme- 
diate redshift to the local value , we also chose to normalize to 
the galaxy B band luminosity. 

The local SN rates measured by collecting d ata of five pho- 
tographic SN searches (Cappell aro et alj [1999: and reference 
therein) are: r Ia = 0.17 + 0.04 and r cc = 0.41 + 0.17 /z^SNu 
at z — 0.01. As it can be seen the SN la rate in SNu appears con- 
stant, within the uncertainties, up to z — 0.3, whereas the SN CC 
rate increases by a factor 2 already at z = 0.21. The different 
evolutionary behaviour of CC and la SN rates implies that their 
ratio increases by a factor of ~ 2 from the local Universe to a 
look-back time of "only" 3 Gyr (z = 0.25): 

- (r cc /r / >. i ) =2.4 ±1.1 

- (r cc /r la \ . 25) =5.6 + 3.5 

Considering that, in this same redshift range, the cosmic SFR 
nearly doubles, the evolution with redshift of the ratio r cc /r Ia 
requires that a significant fraction of SN la progenitors has a 
delay time longer that 2-3 Gyr (cf. Sect. |8j. 

The interpretation on the evolution of the rate in SNu is not 
straightforward, as it reflects both the redshift dependence of 
both the SFH and the B-band luminosity. Indeed, the Z?-band 



luminosity, with contribution from both old and young stars, 
evolves with a different slope in comparison with the on-going 
SFR. 

We acknowledge that the estimate of the SN rate evolution 
with redshift depends on the adopted extinction correction. For 
instance, the ratio between the rate at z=0.21 and the local rates 
for CC SNe varies from 1.6 to 2.8, when no or high extinction 
correction is applied. Anyway, the fact that the CC SN rate in- 
creases faster than the SN la rate appears a robust result. 

7.4. SN rates and galaxy colors 

Integrated colours are valuable indicators of the stellar popula- 
tion and SFR in galaxies and, for galaxies at high redshift, the 
colour information is easier to derive than the morphological 
type. Therefore it is interesting to investigate the dependence of 
the SN rates on the galaxy colours. 

In the local Universe , SN rates as a f unctio n of galaxy 
colours were derived by |Cappellaro et al.1 dl999f) for optical 
bands and by iMannucci et all d2005l) for optical-infrared bands. 
The SN la rate per unit B luminosity appears almost constant 
in galaxies with different U — V color, whereas the CC SN rate 
strongly increases from red to blue galaxies, a trend very similar 
to that of the SFR. This indicates that a fair fraction of SN la pro- 
genitors has long delay times, whereas CC SNe are connected to 
young massive stars. 

On the other hand, when the local SN rates are normalized 
to the K band luminosity, or galaxy mass, they show a rapid in- 
crease with decr easing B - K galaxy color for all the SN types 
(Mannucci et al. 2005). This result indicates that a sizeable frac- 
tion of SNe la has short delay times. A population of SN la pro- 
genitors with short delay time has been propo sed also to explain 
the h igh S N la rate in radi o -loud galaxies dDella Valle et al. I 
120051) . and IMannucci et al. I (120061) suggest that this fraction 
("prompt" component) amounts to ~ 50%. In this respect we 
notice that evolutionary scenarios for SN la progenitors do in- 
deed predict a distribution of the delay times ranging from a few 
tens of Myr to the Hubble time, or more, and that su ch distri- 
butions are typically skewed at the early delay times (Greggio 
Il983t|2005)." Given the diagnostic power of the trend of the SN 
rates with the colour of the parent galaxy we investigate here this 
trend at intermediate redshift. 

In the last decade the extensive study of large samples 
of galaxies at different redshifts has sho wn the existence of 
a bimodal distribu t ion of galaxy c olors dStrateva et al. I l200~lt 
IWeineret all 120051: iBell et all 120041) . In both apparent and rest 
frame color magnitude diagrams, galaxies tend to separate into 
a broad blue sequence and a narrow red sequence. While the red 
sequence is adequately reproduced by passive evolution mod- 
els, the blue sequence hosts galaxies with different rate of on- 
going star formation . Actually, a minor fraction of the red galax- 
ies are also star forming (e.g. edge-on disks or starbusts), their 
colours being red colours because of high dust extinction. The 
fraction of star formi ng galaxies in the red population is 15-20% 
dStrateva et al. 1 1200 lb at low redshift, and i ncreases to 30% at 
z ~ 0.7 dBell et al.ll2004HWeiner et al.ll2005l) . The bimodality is 
not restricted to colour alone, but extends to many other galaxy 
properties as luminosity, mass and environment, which are all 
well correlated with colours. 

We split up our galaxy sample into blue and red sub-samples, 
according to the observed B - V color. We took the rest frame 
color of Sa CWW temp l ate as reference. The local galaxy sam- 
ple of Cappell aro et a U d 19991) was divided in the same way. 
U — V or B — R colours could be more sensitive tracers of the 
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Table 10. SN rates as a function of galaxy color [SNu h 2 ] 



SN type 


r(z = 

red 


0.01) 
blue 


r(z = 
red 


0.25) 
blue 


la 


0.20 ± 0.04 


0.19 ±0.04 


16 +0 13 

-0.18 


23 +013 

u -0.09 


CC 


0.10 ±0.04 


0.99 ±0.15 


57 +0 - 48 

-0.30 


1 50 +0 - 66 


All 


0.30 ± 0.06 


1.18 ± 0.16 


73 +a5 ° 


i 70 +0.67 
±./0_ n4 , 



stellar populations, but [/-band photometry is not available for 
our galaxy sample and /?-band is not available for the local sam- 
ple. 

The SN rates in the blue and red galaxy samples were com- 
puted by distributing the unclassified SN and SNAGN candi- 
dates among type la and CC SNe, based on the observed frac- 
tions of the spectroscopically confirmed SNe in each of the two 
galaxy samples (in the red sample SNe la are 67%, whereas they 
are only 25% in the blue sample). 

The SN rates (in SNu) in the local Universe and at redshift 
Z = 0.25 for the red and blue galaxy subsamples are listed in 
Tab. [TO] 

We found very similar trends in the local Universe and at red- 
shift z = 0.25 both for type la and CC SNe. It appears that while 
the SN la rate is almost costant in galaxies of different colors, 
the CC SN rate alw ays peaks in blue ga laxies. This result is con- 
sistent with that of ISullivan et al. I (120061) who, af ter comparison 
with the measurements of M annucci et al. (2005), found that the 
SN la rate as a function of galaxy colours does not evolve sig- 
nificantly with redshift. 

The rapid increase with redshift of the CC SN rate both in 
blue and red galaxies can be attributed to an increasing propor- 
tion of star forming galaxies in the red sample going from low 
to intermediate redshift. 



8. Discussion 

In this section we discuss the evolution of the SN rates and in- 
vestigate on the link between SF and SN rates. First we ver- 
ify the consistency of our estimate of the redshift evolution of 
SN rates with other measurements from the literature (Sec. 18. U . 
Then we compare the observed evolution of the CC SN rate 
with those expected for different SFHs (Sec. I8.2t . Finally, we 
convolve the SFH o f lHopki ns & Beacom I (2006) with the delay 
time distribution (DTD) for different progenitor scenarios, and 
compare the results to the measurements of the SN la rate evo- 
lution (Sec. 



8.1. Comparison with other estimates 

Published measurements of the SN rates at intermediate and high 
redshift are expressed in units of co-moving volume. To convert 
our rates from SNu to volumetric rates we assumed that the SN 
rates are proportional to the galaxy B band luminosity. If this 
assumption is true, by multiplyng the rates in SNu by the total B 
band luminosity density in a given volume we derived the total 
rates in that volume, even if we did not sample the faint end of 
the galaxy luminosity function. We accounted for the evolution 
of the B band luminosity density with redshift. A compilation of 
recent measurements of the B band luminosity density is plotted 
as function of redshift in Fig [2T| where we also show a linear 
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Fig. 21. Measurements of the galaxy luminosity density at differ- 
ent redshift s. The DEEP 2 and COMBO-17 data are taken from 
Table 2 in iFaber et "ail (120051) . The line represents the linear 
least-square fit in the redshift interval z = - 1 . 



least-square fit to the data in the range z — 0— 1: ]b{z) = (1.03 + 
1.76 xz) 10 8 Z|Mpc~ 3 . 

Multiplying our measurements by the value of js at the av- 
erage redshifts of the la and CC SN samples, the rates per unit 
of co-moving volume result: 

- r ld (z = 0.30) = 0.34+°;J£°;g lo-^yr-'Mpc" 3 

- r cc (z = 0.21) = 1.15!°;™2 lO^yr^Mpc- 3 

where both statistical and systematic errors are indicated. Also 
local rates are converted into volumetric rates: r Ia (z = 0.01) = 
0.18 + 0.05 lO-^yr^Mpc -3 and r cc (z = 0.01) = 0.43 + 0.17 
10" 4 /! 3 yr _1 Mpc -3 . With respect to the rates in SNu's, the vol- 
umetric rates evolve more rapidly with redshift, due to the in- 
crease of the B band luminosity density. We find an increase of a 
factor of ~ 2 at z — 0.3 for SNe la, and a factor of ~ 3 at z — 0.21 
for CC SNe. 

Measurements of la and CC SN rate as function of red- 
shift are shown in Fig. l22l where those originally given in SNu 
(lHardin et al. Il2000t (Blanc et al. 112001: ICappellaro etai7ll2005l) 
were converted into measurements per unit volume as above. 

As it can be seen, the few measurements of the CC SN rate 
appear to be fully consistent, while those of the SN la rate show a 
significant dispersion which increases with redsh ift, in particular 
in the rang e 0.5 < z < 0.7 wher e the values of Barri s & Tonry I 
(2006 1) andlDahlen et all d2004l) are 2-3 times higher that those 
of iPain et al. I (120021) and lNeill et al.l (120071) . Our estimate of the 
SN la rate is consistent with all other measurements in the red- 
shift range we explored; our result does not help to discriminate 
bet ween the steep t rend s uggested by the Barris & Tonry (2006) 
and Dahlen et al. 1 (120041) measurements and the slow evolution 
indicated by fhe lNeill et al.l (120071) measurement. The robust in- 
dication from the current data appears that the SN la rate per 
unit volume at redshift 0.3 is a factor of ~ 2 higher that in the 
local Universe, while in the same redshift range the CC SN rate 
increases by a factor of ~ 5. 
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Fig. 22. Observed SN rates as function of redshift from different authors as indicated in the legend. The black (gray) symbols 
indicate SN la (SN CC) rate measurements. The shaded area represents the 1 cr confidence level of our rate evolution estimate as 
deduced from the MLE fit. 



8.2. Comparison with the predicted evolution of the CC SN 
rate 



The stellar evolution theory predicts that all stars more massive 
than 8-10 M complete the eso-energetic nuclear burnings, up to 
the development of an iron core that cannot be supported by any 
further nuclear fusion reactions or by electron degenerate pres- 
sure. The subsequent collapse of the iron core results into the 
formation of a compact object, a neutron star or a black hole, ac- 
companied by the high-velocity ejection of a large fraction of the 
progenitor mass. Typell SNe originate from the core collapse of 
stars that, at the time of explosion, still retain their H envelopes, 
whereas the progenitors of type Ib/Ic SNe are thought to be mas- 
sive s tars which have lost their H (and He) envelope (Heg er et al.l 
120031) . Given the short lifetime of their progenitors (< 30 Myr), 
there is a simple, direct relation between the CC SN and the cur- 
rent SF rate: 



where if/(z) is the SFR and K cc is the number of CC SN 
progenitors from a 1 M Q stellar population: 



££c <p{m)dm 
J™ u m(f>(m)dm 



(8) 



r CC {z) 



(7) 



where 0(m) is the IMF, mi - mu is the total stellar mass range, 
and mf c - m^ c is the mass range of CC SN progenitors. 

In principle, if accurate measurements of the CC SN and SF 
rates are available, it is possible to probe the possible evolution 
with redshift of either the CC SN progenitor scenarios or the 
IMF by determining the value of K cc . Here however we assume 
that K cc does not evolve significantly in the redshift range of 
interest and compare the observed evolution of the CC SN rate 
with that predicted by the SFH. 

The estimate of the cosmic SFH is based on different SF in- 
dicators, depending on redshift range, and requires a suitable 
parameterization and accurate normalization. Since there is a 
large scatter between the measurements obtained by different 
SFR indicators (a factor 2-3 at z = 0.3 - 0.5 and increasing 
with redshift) it is hard to o btain a consistent picture of the SFH 
dHopkins & Be acom 2006). Indeed observations made at differ- 
ent wavelengths, from X-rays to radio, sample different facets 
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Fig. 23. Comparison between the SN CC and SF rate evolution. 
Symbols are as in Fig. [22] with additional open symbols (mea- 
surements not corrected for extinction) and filled black symbols 
(estimates for the high extinction correction). Lines are selected 
SFR evolutions from the literature. All SFHs have been scaled 
to the SalpeterA IMF. 

of the SF activity and are sensitive to different time scales over 
which the SFR is averaged. Thus, different assumptions are re- 
quired to convert the observed luminosities at the various wave- 
lengths to the SFR, and different systematic uncertainties affect 
the SFR estimates. In particular the significant difference be- 
tween the SFR inferred from the UV and Ha luminosity and 
that inferred from the far-IR (FIR) luminosity may be related to 
the effect of dust extinction (expecially for the UV) and/or to the 
contribution to the light from old stars and AGNs (for the FIR). 

To illustrate this point, we select three representative pre- 
scriptions for the SFH, namely: 

- the piecewise linea r fit of selected SFR measu rements in the 
range < z < 6 by Hop kins & Beacom I d2006l) . 

- the fit to the SFR measurem ents from the Ha emission 
line, by Hippel ein et al. Id2003l) . with an exponential increase 
from z = to z = 1 .2, 

- the prescriptions bv lHernquist & Springei1 (l2003l) . i.e. a dou- 
ble exponential function that peaks at redshift z ~ 5.5, ob- 
tained from an analytical model and hydro-dynamic simula- 
tions. 

All these SFHs were converted to the same IMF, a modi- 
fied Salpeter IMF (SalA) with a turnover below 0.5 M and 
defined in the mass range m i - 0.1M o to mu = 120 M 
dBaldrv & Glazebrook 1 [2003). We assumed a mass range of 
8 - 5OM for CC SN progenitors, which gives a scale factor 
K cc = 0.009. 

The measurements of the CC SN rate per unit volume and the 
predicted evolutionary behaviours are shown in Fig [23] The ob- 
servations con firm the steep increase with reds hift expected by 
the SF H from iHopkins & Beacoml d2006l) and iHippelein et al. I 
(2003;). For a look-back time of 3 Gyr (z = 0.25) both the 



SFR and the CC SN rate increase by a factor of ~ 3 com- 
pa red with the local values. A fl at evolution, as that proposed 
by iHernquist & Springel I d2003l) . appears inconsistent with the 
observed CC SN rates in the overall range of redshift. 

With the ad opted K cc , the level of t he CC SN rate predicted 
by the SFH of Hip pelein et al. 1 (120031) and, in general, by the 
U V and Ha based SFHs fits w ell the data. Instead, for the SFH 
of Hopkin s~& BeaconTI d2006l) and in general the SFHs inferred 
through FIR luminosity, the predicted CC SN rate is higher than 
observed over the entire redshif t range (see also Dahlen et al. I 
I2004t IHopkins & Beacom1l2006t iMannucci et al. |[2007l) . 

If we correct the observed CC SN rate according to the high 
extinction scenario, we obtain an a cceptable agreement between 
the data and the predictions with Hopkin s~& BeaconTI ([2006) 
SFH, as shown in Fig. [23] However, this correction would re- 
quire an extremely high dust content in galaxies which is not fa- 
vored by present measurements. Indeed, M annucci et al. I {2007) 
derived an estimate of the fraction of SNe which are likely to be 
missed in optical SN searches because they occur in the nucleus 
of starburst galaxies or, in general, in regions of very high ex- 
tinction and found that, at the average redshift of our search, the 
fraction of missing CC SNe is only ~ 10%, far too small to fill 
the gap between observed and predicted rates. 

Alternatively we may consider the possibility of a narrower 
range for the CC SN progenitor masses: in particular, a lower 
limit of 10 - 12 M would bring the observed CC SN rates in 
agreement with those predicted from FIR based SFHs. In this 
respect we notice that, from a theoretical point of view, there is 
the possibility that a fraction of stars between 7-8 M and 10- 
12 M avoids t he collapse of the core and ends up as ON eMg 
White Dwarfs dRitossa et al .11 19981: iPoelarends et alj|2007l) . On 
the other hand, estimates of the progenitor mass from the de- 
tection in pre-explosion (HST) images has been p ossible for a 
few SNe IIP (e.g. SN 2 003gd dHendrv et al.ll2005h . SN 2005cs 
dPastorello et al. Il2006l) ): their absolute ma gnitudes and colours 
seem indicate a moderate mass (8- 12 Mn )dVan Dvk et al.l2003h 
ISmartt et al. 112004 iMaund et alJ2005tlLi et al. Il2006l) ~ 

Given these controversies, we conclude that in order to con- 
strain the mass range of CC SN progenitors it is necessary to 
reduce the uncertainties in the cosmic SFH. In addition it is im- 
portant to apply a consistent dust extinction correction both to 
SFH and to CC SN rate. 



8.3. Comparison with the predicted evolution of SN la rate 

According to the standard scenario SNe la originate from the 
thermonuclear explosion of a Carbon and Oxygen White Dwarf 
(C-O WD) in a binary system. In the first phase of the evolution 
the primary component, a star less massive than 8M , evolves 
into a C-O WD. When the secondary expands and fills its Roche 
Lobe, two different paths are possible, depending on whether a 
common envelope forms around the two stars (double degener- 
ate scenario, DD) or not (single degenerate scenario, SD). In the 
SD scenario, the WD remains confined within its Roche Lobe, 
grows in mass until it reaches the Chandrasekhar limit and ex- 
plodes, while in the DD scenario the binary system evolves into 
a close double WD system, that merges after orbital shrinking 
due to the emission of gravitational wave radiation. 

In both scenarios two basic ingredients are required to model 
the evolution of the SN la rate: the fraction of the binary sys- 
tems that end up in a SN la, and the distribution of the time 
elapsed from star formation to explosion (delay-time). In the 
SD scenario, the delay time is the evolutionary lifetime of the 
secondary; in the DD, the gravitational radiation timescale has 
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to be added. In both cases, the distribution of the de lay times 
depen ds on the distribution of the binary parameters dGreggio I 
120051) . In principle, the SD and DD scenarios correspond to dif- 
ferent realization probabilities and different shape of the delay 
time distribution (DTD) functions, hence rather different evolu- 
tionary behavior of the SN la rate. As mentioned previously, the 
observations indicate that the distribution of the delay times of 
SN la progenitors is rather wide. 

The cosmic SFH is the other critical ingre dient that modu- 
lates the evolution of the SN la rate. Following lGreggiol f2005) 
the SN la rate is given by: 



r Ia (f) = k 



r 

aA J 71 



•min(t,T x ) 



(9) 



where k a is the number of stars per unit mass of the stellar 
generation, A Ia is the realization probability of the SN la scenario 
(the number fraction of stars from each stellar generation that 
end up as SN la), / Ia (r) is the distribution function of the delay 
times and if/(t - r) is the star formation rate at the epoch t — t. The 
integration is extended over all values of the delay time t in the 
range t, and min(t, r x ), with r, and t x being the minimum and 
maximum possible delay times for a given progenitor scenario. 
Here we assumed that both k a and A Ia do not vary with cosmic 
time. 

A detailed analysis of the predicted evolution of the 
SN la r ate for differ e nt SFHs and DTDs is p r esente d else- 
where dForsteret alJ 120061: iBlanc & Greggiol 120071) . Here 
we consider only one SF H, the piecewise interpolation of 
Hop kins & Beacoml d2006l) . since this is conveniently defined 
also at high redshift, and limit our analysis to few DTDs repre- 
sentative of different appr oaches to model the SN la rate evo- 
lution: three models from Greggio (2005 ), and two differen t 
parametrizations (Mann ucci et al. 1 120061 : IStrolger et al. I [2004) 
designed to address some specific observational constraints, re- 
gardless the correspondence to a specific progenitor scenario. 

Specifically, among the Greggio (2005) models we select 
one SD and two DD-models, one of the "close" and the other of 
the "wide" variety, the latter being an example of a relatively fiat 
DTD. This choice is meant to represent the full range of plau- 
sible DTDs. The minimum delay time for the SD model is the 
nuclear lifetime of the most massive secondary stars in the SN 
la progenitor's system, i.e. 8M . In principle, for the DD model 
the minimum delay time could be appreciably larger than this 
because of the additional gravitational waves radiation delay. In 
practice, also for the DDs the minimum delay is of a few 10 7 yrs. 
The maximum delay time is quite sensitive to th e model for the 
SN la progenitors. The reader should refer to iGreggiol 12005) 
for a more detailed description of these models. 

Th e DTD parametrization proposed by iMannucci et al. I 
(2006) is the sum of two distinct functions: a Gaussian centered 
at 5xl0 7 yr, representative of a "prompt" progenitor population 
which traces the more recent SFR, and an exponentially declin- 
ing function with characteristic time of 3 Gyr, a "tardy" progen- 
itor population proportional to the total stellar mass. This DTD 
was introduced to explain, at the same time, the dependence of 
the SN la rate per unit mass on the galaxy morphological type, 
the cosmic evolution of the SN la rate and, in particular, the high 
SN la rate observed in radio-loud galaxies. 

Finally, IStrolger et aTl d2004 showed that the best fit of 
the apparent decline of SN la rate at z > 1 is achieved for 
a DTD with a Gaussian distribution centered at about 3 Gyr. 
Nevertheless this DTD fails to reproduce the dependence of 
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Fig. 25. SN la rate me asurements fitted with diffe rent DTD func- 
tions and the SFH by Hopki ns~& Beacoml d2006l) . Symbols for 
measurements are as in Fig. [22] 



the SN r ate on galaxy colours which is observed in the local 
Universe (Mannucc i et al. Il2006l) . 

The selected DTD functions are plotted in Fig. [24] while the 
predicted evolutionary behaviours of the SN la rate are com- 
pared with all published measurements in Fig. [25] In all cases, 
the value of £ (! A Ia was fixed to match the value of the local rate; 
depending on the model it ranges between 3.4-7.6 xl0~ 4 . This 
normalization implies that, for the adopted SalA IMF, and as- 
suming a mass range for the progenitors of 3 - 8 M Q , the proba- 
bility that a star with suitable mass becomes a SN la, is ~ 0.01- 
0.03. 

The models obtained with the different DTDs are all con- 
sistent with the observations with the exception of the "wide" 
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DD model, whose redshift evolution is definitely too flat. On 
the other hand none of the DTD functions, with the adopted 
SFH, is able to reproduce at the same time the very rapid in- 
crease from redshift t o 0.5 suggested by some measurements 
dBarris & Tonrv 1 120061: IDahlen et : al. 1 12004 and the decline at 
redshift > 1 (IDahlen et al. 20041) . We note that a new mea- 
surement of Poznanski et al. (2007) suggests that the SN la rate 
decline at high red shift may be not as steep as estimated by 
IDahlen et al. Id2004l) . 

Given the current uncertainties of both SN la rate and SFH it 
is difficult to discriminate between the different DTD functions 
and hence between the different SN la progenitor models. To 
improve on this point, more measurements of SN la rate at high 
redshifts are required to better trace the rate evolution. At the 
same time measurements in star forming and in passive evolving 
galaxies in a wide redshift range can provide important evidence 
about the SN la progenitor models. In addition, it is essential to 
estimate the cosmic SFH more accurately because the position 
of the peak of the SFH w as found to be the c rucial parameter for 
the recovered delay time dForsteret al.ll2006h . 



9. Summary 

In this paper we have presented our SN search (STRESS) carried 
out with ESO telescopes and aimed at measuring the rate for 
different SN types at intermediate redshift. 

Our approach, that consists in counting the SNe observed in 
a well defined galaxy sample, allows us to investigate the de- 
pendence of the rates on galaxy colours, and to perform a direct 
comparison with measurements in the local Universe, which are 
obtained with the same approach. 

For the selection and the characterization of the galaxy sam- 
ple we used multi-band observations and the photometric red- 
shift technique. After 16 observing runs, 12 in V, 4 in R band, 
we discovered a few hundred variable sources from which we 
selected 86 SN candidates: 25 spectroscopically confirmed SNe 
(9 SNe la and 16 CC SNe), 33 SN candidates and 28 SNAGN 
candidates. 

The control time of our galaxy sample has been evaluated 
through a detailed analysis of the SN detection efficiency of our 
search, with a number of MonteCarlo simulations. In addition, 
we have improved the handling of the correction for dust extinc- 
tion by elaborating realistic modelling. The rates are presented 
for two scenarios, one stardard and one which, maximizing the 
effect, yield reasonable upper limits. 

Since our SN sample spans a wide redshift range (0.05 < 
Z < 0.6) we have obtained an observational constraint on the 
redshift evolution of the rates, described as a power law with 
two parameters: r(z), the rate at the mean redshift of the search 
z, and a, an evolution index. For SNIa and CC SNe we get: 

SN la: r(z = 0.3) = 022 + ^_° Q \ 6 4 h 2 10 SNu, a = 4.4^3.3 
CC SN: r(z = 0.2) = 0.82^;™° h 2 JQ SNu, a = 7.5^-8+3.6 

where both statistical and systematic uncertainties are re- 
ported. 

Our results indicate that, compared to the local value, the CC 
SN rate per unit B band luminosity increases by a factor of ~ 2 
already at z ~ 0.2, whereas the SN la rate is almost constant up to 
redshift z ~ 0.3. The dependence of the SN rates per unit B band 
luminosity on the galaxy colour is the same as observed in the 
local Universe: the SN la rate seems to be almost constant from 
red to blue galaxies, whereas the CC SN rate seems to peak in the 
blue galaxies. Therefore, on the one hand, the observed evolution 



with redshift of the ratio r cc jr la requires a significant fraction of 
SN la progenitors with long delay times; on the other hand, the 
observed trend with galaxy colour requires a short delay time for 
a fraction of SN la progenitors. 

Systematic uncertainties in our SN rate measurements have 
been investigated by considering the influence of the different 
sources. In particular we have analysed the important role played 
by the AGN contamination of SN candidate sample and the host 
galaxy extinction. 

After accounting for the evolution with redshift of the blue 
luminosity density, we have compared our estimates with all 
other measurements at intermediate redshift available in the lit- 
erature and found that they are consistent, within the relatively 
large errors. 

Finally we have exploited the link between SFH and SN rates 
to predict the evolutionary behaviour of the SN rates and com- 
pare it with the path indicated by observations. 

The predicted evolution of the CC SN rate with redshift 
has been computed assuming three representative SFHs and the 
mass range 8-50 M Q for CC SN progenitors. The comparison 
with the observed evolution confirms the steep increase with red- 
shift indicated by most recent SFH estimates. Specifically, we 
found a good agreement with the predictions from t he SFH in- 
ferred through Ha luminosity dHippelein et al.] [2003), while the 
SFHs from FIR luminosity dHopkins & Beacom 1 2006) overesti- 
mate the CC SN rate, unless a higher extinction correction and/or 
a narrower range for the progenitor masses is adopted. This re- 
sult illustrates how interesting clues can be obtained by compar- 
ing the SN CC rate to other SFR tracers, in the same galaxy sam- 
ple, to verify the reliability of the techniques used to derive SFR 
estimates and the adequacy of the dust extinction corrections. 

The cosmic evolution of the SN la rate ha s been estimated 
by convolving the SFH of Hopki ns"& Beacoml d2006l) with var- 
ious formulations of DTD related to different SN la progenitor 
models. All DTDs appear to predict a SN la rate evolution con- 
sistent with the observations, with the exception of the "wide" 
DD model, which appears too flat. At the same time none of 
the explored DTD functions, at least with the adopted SFH and 
k a A l,i factor are able to reproduce simultaneously both the rapid 
increase from redshift to 0.5 and the decline at redshift > 1 
suggested by some measurements. 

With the current data of the rate evolution, it is difficult to 
discriminate between different DTDs and then between different 
SN la progenitor models. Measurements of the SN la rate in 
star forming and in passive evolving galaxies in a wide range 
of redshifts can provide more significant evidence about SN la 
progenitors. 

An extensive survey to search all SN types in a well char- 
acterized sample of galaxies is of the highest priority to probe 
the link between the SN rates and the SFR in different contexts. 
This will yield important constraints on SN progenitor scenarios 
from a detailed analysis of SN rates as a function of the various 
properties of galaxies. 
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Table 2. Observing log of the SN search fields. 



run 


Field 


Filter 


Nexp 


Texp 


Seeing 


1999-02-23 


10Z2 


B 


3 


1200 


1.12 




13Z3 


B 


3 


1200 


1.10 


1999-03-10 


10Z2 


B 


3 


1200 


1.14 




13Z3 


B 


3 


1200 


1.15 




10Z2 


V 


3 


1200 


1.25 




13Z3 


V 


3 


1200 


1.15 


1999-03-19 


10Z2 


B 


3 


1200 


1.04 




13Z3 


B 


3 


1200 


0.96 




10Z2 


V 


3 


1200 


0.86 




13Z3 


V 


3 


1200 


0.84 


1999-05-08 


10Z2 


B 


3 


1200 


1.34 




13Z3 


B 


3 


1200 


1.4 




Field2 


B 


3 


1200 


1.91 




Field2 


R 


3 


1200 


1.02 




10Z2 


V 


3 


1200 


1.37 




Field2 


V 


3 


1200 


1.62 


1999-05-17 


10Z2 


B 


3 


1200 


0.81 




13Z3 


B 


3 


1200 


1.74 




Field2 


B 


3 


1200 


0.74 




10Z2 


R 


3 


1200 


0.96 




13Z3 


R 


3 


1200 


0.91 




10Z2 


V 


3 


1200 


0.90 




13Z3 


V 


3 


1200 


0.98 




Field2 


V 


3 


1200 


0.81 


1999-08-03 


13Z3 


B 


3 


1200 


1.43 




22Z3 


B 


3 


1200 


1.17 




Field2 


B 


3 


1200 


1.12 




J1888 


B 


3 


720 


1.22 




J1888 


R 


3 


720 


1.31 




Field2 


V 


3 


1200 


1.8 




J1888 


V 


3 


720 


1.41 


1999-08-14 


22Z3 


B 


2 


1200 


1.58 




EisAl 


B 


3 


1200 


1.67 




Field2 


B 


3 


1200 


2.15 


1999-08-31 


22Z3 


B 


3 


1200 


1.43 




Field2 


B 


3 


1200 


2.44 


1999-09-13 


J1888 


B 


5 


720 


1.62 




J1888 


R 


5 


720 


1.24 




22Z3 


V 


3 


1200 


1.46 




EisAl 


V 


3 


1200 


1.32 




EisB 


V 


3 


1200 


1.27 




Field2 


V 


3 


1200 


1.67 




J1888 


V 


5 


720 


1.41 


1999-10-30 


AXAF 


B 


3 


900 


3.00 




EisAl 


B 


3 


900 


1.32 




EisB 


B 


3 


900 


2.02 




EisAl 


V 


3 


900 


1.50 




EisB 


V 


3 


900 


1.52 


1999-11-09 


03Z3 


B 


3 


900 


0.98 




AXAF 


B 


3 


900 


1.13 




EisAl 


B 


3 


900 


1.12 




J1888 


B 


3 


900 


1.05 




03Z3 


V 


3 


900 


1.07 




AXAF 


V 


3 


900 


1.08 




EisAl 


V 


3 


900 


0.84 




EisB 


V 


3 


900 


0.76 




J1888 


V 


3 


900 


0.81 


1999-12-02 


AXAF 


B 


3 


900 


1.03 




AXAF 


V 


3 


900 


1.12 




EisB 


V 


3 


900 


1.27 
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Table 2. continued. 



run 



1999-12-10 



1999-12-28 



2000-01-08 



2000-11-16 



2000-12-17 



2001-04-18 



2001-04-19 



2001-11-11 



2001-11-12 



2001-11-18 



Field 


Filter 


Nexp 


Texp 


Seeinj 


EisC 


V 


3 


900 


1.15 


EisD 


V 


3 


900 


1.2 


J1888 


V 


3 


900 


1.06 


03Z3 


B 


3 


900 


1.6 


03Z3 


V 


3 


900 


1.62 


10Z2 


V 


3 


900 


1.28 


EisAl 


V 


3 


900 


1.75 


EisC 


V 


3 


900 


1.31 


J1888 


V 


3 


900 


1.65 


EisC 


B 


3 


900 


1.18 


03Z3 


V 


3 


1200 


2.18 


AXAF 


V 


3 


900 


1.15 


EisC 


V 


3 


900 


0.77 


EisD 


V 


3 


900 


0.96 


J1888 


V 


3 


900 


1.10 


EisC 


B 


3 


900 


1.06 


EisD 


B 


3 


900 


0.91 


EisC 


V 


3 


900 


0.85 


EisD 


V 


3 


900 


0.90 


10Z2 


V 


3 


900 


0.92 


03Z3 


V 


3 


900 


1.14 


AXAF 


V 


3 


600 


0.92 


EF0427 


V 


3 


900 


0.74 


EisAl 


V 


3 


900 


0.92 


EisB 


V 


3 


900 


0.86 


EisC 


V 


3 


900 


0.64 


EisD 


V 


3 


900 


0.83 


J1888 


V 


3 


900 


1.05 


03Z3 


V 


3 


600 


1.14 


AXAF 


V 


3 


600 


0.93 


EF0427 


V 


3 


900 


0.93 


EisAl 


V 


3 


900 


1.12 


EisB 


V 


3 


900 


1.26 


EisC 


V 


3 


900 


0.85 


J1888 


V 


3 


660 


1.02 


10Z2 


V 


3 


900 


0.98 


13Z3 


V 


3 


900 


0.85 


22Z3 


V 


3 


900 


1.41 


EisC 


V 


3 


900 


0.85 


EisD 


V 


3 


900 


0.86 


Field2 


V 


3 


900 


0.84 


white23 


V 


3 


900 


0.82 


white27 


R 


3 


1200 


0.71 


whitehz2 


R 


3 


1200 


0.67 


white27 


V 


3 


900 


0.73 


white3 1 


V 


3 


900 


1.01 


whitehz2 


V 


3 


900 


0.79 


03Z3 


V 


3 


900 


0.90 


22Z3 


V 


3 


900 


0.84 


EF0427 


V 


3 


900 


0.75 


EisAl 


V 


3 


900 


0.82 


EisB 


V 


3 


900 


0.78 


EisC 


V 


3 


900 


0.91 


J1888 


V 


3 


900 


0.75 


22Z3 


R 


3 


900 


1.31 


AXAF 


R 


3 


900 


0.92 


EisAl 


R 


3 


900 


1.55 


EisB 


R 


3 


900 


1.07 


J1888 


R 


3 


900 


1.24 


AXAF 


V 


3 


900 


1.03 


03Z3 


V 


3 


900 


0.90 
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Table 2. continued. 



run 



Field 


Filter 


Nexp 


Texp 


Seein, 


22Z3 


V 


3 


900 


1.10 


AXAF 


V 


3 


900 


0.84 


EF0427 


V 


3 


900 


0.93 


EisAl 


V 


3 


900 


1.20 


EisB 


V 


3 


900 


0.83 


EisC 


V 


3 


900 


0.65 


J1888 


V 


3 


900 


1.02 


03Z3 


V 


3 


900 


0.94 


AXAF 


V 


3 


900 


0.95 


EF0427 


V 


3 


900 


0.95 


EisB 


V 


3 


900 


0.88 


EisC 


V 


3 


900 


0.95 


J1888 


V 


3 


900 


0.94 


22Z3 


R 


1 


900 


2.00 


AXAF 


R 


3 


900 


0.82 


EF0427 


R 


3 


900 


1.1 


EisB 


R 


3 


900 


0.89 


EisC 


R 


3 


600 


1.05 


J1888 


R 


3 


900 


0.84 


22Z3 


V 


1 


900 


1.90 


EisD 


V 


3 


900 


0.92 


EisC 


V 


3 


900 


1.32 


EisD 


V 


3 


900 


0.92 


10Z2 


V 


3 


900 


0.91 


13Z3 


V 


3 


900 


0.78 


white27 


V 


3 


900 


0.80 


white3 1 


V 


3 


900 


0.85 


Field2 


V 


3 


900 


1.30 


EisD 


R 


3 


900 


0.85 


10Z2 


R 


3 


900 


0.88 


13Z3 


R 


3 


900 


1.1 


white3 1 


R 


3 


900 


0.87 


Field2 


R 


3 


900 


1.31 


white23 


V 


3 


900 


0.96 


whitehz2 


V 


3 


900 


0.91 


EF0427 


R 


3 


900 


0.85 


EisC 


R 


3 


900 


0.91 


EisD 


R 


3 


900 


0.93 


white27 


R 


3 


900 


0.82 


white3 1 


R 


3 


900 


0.94 


white3 1 


R 


3 


900 


0.91 


whitehz2 


R 


3 


900 


0.83 


10Z2 


R 


3 


900 


0.87 


EF0427 


I 


3 


900 


1.05 


EisC 


I 


3 


900 


0.80 


EisD 


I 


3 


900 


1.03 


10Z2 


I 


3 


900 


0.73 


whitehz2 


I 


3 


900 


0.76 


white27 


I 


3 


900 


0.67 


white3 1 


I 


3 


900 


0.75 


13Z3 


R 


3 


900 


0.78 


EisD 


R 


3 


900 


1.11 


22Z3 


R 


3 


900 


0.86 


white3 1 


R 


3 


900 


1.18 


whitehz2 


R 


3 


900 


0.78 


white27 


R 


3 


900 


0.81 


EisD 


B 


3 


900 


1.57 


white3 1 


B 


3 


900 


1.38 


whitehz2 


B 


3 


900 


1.38 


white27 


B 


3 


900 


1.55 


white27 


R 


3 


900 


1.01 



2001-12-08 



2001-12-09 



2002-04-07 



2002-04-08 



2004-02-18 



2004-02-19 



2004-05-14 



2004-05-15 



2004-07-14 



Botticella et al.: STRESS, Online Material p 5 



Table 2. continued. 



run Field 


Filter 


Nexp 


Texp 


Seeing 


white3 1 


R 


3 


900 


1.12 


22Z3 


R 


3 


900 


0.94 


whitehz2 


R 


3 


900 


0.86 


EisAl 


R 


3 


900 


0.90 


EisB 


R 


3 


900 


1.02 


2005-05-13 white23 


R 


3 


900 


1.27 


wmtez / 


K 


-2 
J 


yuu 


1 A O 


white3 1 


R 


3 


900 


1.55 


whitehz2 


R 


3 


900 


1.18 


13Z3 


R 


3 


900 


1.29 


EisD 


R 


3 


900 


1.56 


10Z2 


R 


3 


900 


1.39 


Field2 


R 


3 


900 


1.32 
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Fig. 1. K correction as a function of SN phase from the observed 
V band to the rest frame B band for redshift 0. 1 . Top panel shows 
the K(BV) for a type la SN, mid panel for a type Ib/c SN and 
bottom panel for a type II SN. In all cases the best fit has been 
obtained neglecting peculiar SNe as SN 1991bg for type la SN 
and SN 1987 A for type II SN. 
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Fig. 15. K(BV) correction for redshift 0.6. 
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Fig. 16. Same as Fig. 1 for the K correction from the observed R 
band to the rest frame B band for redshift 0. 1 . 
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Fig. 14. K(BV) correction for redshift 0.3. 
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Fig. 17. K(BR) correction for redshift 0.3. 
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Fig. 18. K(BR) correction for redshift 0.6. 



